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FOREWORD 


A series of experiments on outdoor propagation of high intensity 
sound is described in this report. Begun in 1975, the research was 
carried out ln_two phases. Phase I, on intense acoustic tones, was 
completed in 1976. Phase II, on intense noise, was finished in 1977. 

The report is primarily devoted to Phase II, but a review of Phase I is 
included. 

Cosponsors of the research along with NASA were Air Force Office 
of Scientific Research, National Oceanic and Atmospheric Administration, 
and Office of Naval Research. 
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SUMMARY 


A series of experiments have been carried out on the propagation of 
finite-amplitude noise outdoors. The purpose of the study was to determine 
the extent to which nonlinear effects influence the propagation of noise 
from a controlled source in the outdoor environment. A secondary purpose 
was to develop, if possible, a theoretical model for the noise propagation. 

A ground-mounted electroacoustic source transmitted broadband, octave 
band, or 1/3 octave band noise in the frequency range 2-10 kHz. The 
source level (overall sound pressure level) of the noise was in the range 
121-145 dB re 20 yPa at 1 m. The propagation path was vertical and parallel 
to an 85 m tower, whose elevator carried the traveling microphone. The 
maximum propagation distance was about 80 m. The experiments were done at 
night during the months June through September 1977. The meteorological 
conditions (ground level) were as follows: temperature range 23-31°C, 

relative humidity range 55-90%, and wind speed range 0-24 km/h. 

The measurements were compared with theoretical predictions based on 
linear theory. Spherical spreading, atmospheric attenuation, and, as 
appropriate, source diffraction were accounted for in these predictions. 
Predictions based on nonlinear theory were also attempted, but the partic- 
ular model was based on an assumption about the noise distortion in the 
transmitter nearfield. As it turned out, this assumption was not justi- 
fied in our experiments. 

Use of an amplitude-frequency scaling law made it possible to compare 
the noise from our experiments with noise from a KC-135A aircraft. The 
importance of nonlinear effects in actual jet noise could thus be estimated. 



The primary conclusions are as follows: 

1. A strong generation of high frequency noise caused by nonlinear 
effects was found in all the high intensity noise experiments. A 
very limited amount of low frequency noise was also generated. The 
intense, middle part of the spectrum deviated little if any from 
expectations based on linear theory. These observations indicate 
that although shocks formed in the noise waveform, distortion did not 
reach the stage at which shock merging was important. 

2. The spectral distortion occurred in both the transmitter nearfield 
and farfield. Moreover, the distortion in the nearfield was over and 
above the spectral changes associated with diffraction. 

3. At no measurement point was small-signal behavior established 
for the high frequency noise. Theoretical calculations for tone 
signals support the proposition that the nonlinearly generated high 
frequency noise never achieves a farfield where small-signal behavior 
is established. 

4. Comparison of the measured spectra with predictions from a 
model based on nonlinear theory showed poor agreement for the high 
frequency noise. The failure of our predictions was not due to an 
error in the nonlinear theory per se, but rather due to an inadequate 
description of the source noise waveform. In the future, use of a 
directly recorded input waveform should be tried. 

5. Comparison of our scaled experimental measurements with actual 
jet spectra show that the spectrum levels encountered in our ex- 
periments are well within the jet noise range. Indeed, the noise 
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measured from a KC-135A jet is roughly 10 dB higher in spectrum 
level than pur scaled noise. One therefore concludes that nonlinear 
effects are probably common in jet noise. 

The report also contains some information about outdoor propagation 
of finite- amplitude tones. Tone experiments done as a forerunner to those 
on noise are reviewed. Additional data on tones taken during the course 
of the noise study are also described. Theoretical analyses of propa- 
gation of very weak and very strong tones are presented. 
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CHAPTER 1 


INTRODUCTION 

Experiments done outdoors on the propagation of Intense acoustic noise 

are the subject of this report. These experiments constitute Phase II of a 

research program begun in June 1975 to study outdoor propagation of finite- 

amplitude acoustics waves, that is, sound so intense that nonlinear 

effects are important. The ultimate application of the research is to jet 

aircraft noise. In particular, the goal has been to determine whether 

aircraft noise is affected by nonlinear propagation distortion. The 

propagation experiments were carried out along an 85 m vertical path 

parallel to a radio tower; the source was on the ground. The first phase 

of the program. Phase I, was done with intense tones as a means of preparing 

for the noise experiments (Phase II) , which were expected to be more 

difficult. Important nonlinear effects — harmonic distortion, shock formation, 

extra attenuation, and an interaction between diffraction and distortion — 

1 2 

were observed in Phase I. * A review is given in Chapter 2. Because of 
the success of the tone experiments, it was decided to proceed with Phase II. 
The results are reported here. 

Previous Investigations of finite-amplitude noise have been largely 

3 _g 9-12 

limited to plane waves in tubes or to strictly theoretical calculations. 

3-5 

Pernet and Payne used a plane wave tube to measure the spectral distortion 
of 1/3 octave and octave bands of cw noise (and, in one experiment, a 2.3 
octave band) . The sound pressure level fePL) at the source was as high as 
140 dB (re 20 yPa) . They observed second and higher harmonic bands of 
noise that grew and decayed with distance according to their theoretical 
predictions. Later, through both theory and experiment, they showed that 
the second harmonic band has an amplitude distribution quite different from 
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8 6 7 

that of the fundamental noise band* Pestorlus ’ 

also made noise propagation measurements In a plane wave tube but used 

generally wider bandwidths and much higher SPL r s, up to 160 dB. In the 

main his noise was pulsed, not -cw. Although we shall presently review 

his results In detail, it may be noted here that because Pestorlus f s 

intensity was greater than Pernet and Payne f s, the low frequency end 

of the spectrum grew substantially as well as the high frequency end. 

Viewed In the time domain, the propagating noise pulse was observed to 

distort into a random sawtooth wave. A computer algorithm based on weak- 

shock theory (with modifications to account for tube wall attenuation and 

dispersion) was used to obtain theoretical predictions. The predictions 

agreed well with the experimental measurements. Later Pestorius's algorithm 

13 

was extended to cover spherical and other nonplanar waves. An attempt 

was then made to use the generalized algorithm to predict the nonlinear 

14 

distortion of actual jet noise. A sample of noise recorded close to a 

jet engine was used as the input to the computer program, and the subsequent 

distortion of the noise signal with distance was calculated. Although 

the calculation did show the importance of nonlinear effects, the-study 

was somewhat artificial because many factors affecting the propagation of 

jet noise in the field — atmospheric attenuation and meteorological effects , 

ground reflection and absorption, and effects due to the spatial extent 

and moving source nature of the jet — were not Included In the algorithm. 

In the meantime, the Russian investigators have been very active on the 
9—12 

theoretical front. They have frequently employed methods developed in 

nonlinear statistical optics. A review of their work on several topics in 
nonlinear distortion of random waves is given in Chap. 10 of Ref. 10. In 
general it is assumed that the initial spectral distribution of the noise 
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and its statistics are known. For several reasons it has been difficult to 
apply the Russian results on broadband noise to our experiments. 

The previous research is sufficient to permit a qualitative description 
of the effect of nonlinearity on noise propagation. The propagation speed 
of a finite^amplitude wave is 

fr - c 0 + b» , (i-i) 

where c q is the small-signal sound velocity, 3 is the coefficient of 
nonlinearity [for gases 3 = (y +l)/2, where y is the ratio of specific heats], 
and u is the particle velocity. Because of its dependence on u, the propa- 
gation speed varies from point to point on the wave. The wave therefore 
distorts as it travels. The compression phases of the wave steepen, and 
the slopes of the expansion phases become more gentle (see, for example. 

Fig. 7-1). Unless the smoothing effect of ordinary absorption is 
stronger than the steepening effect, shoeks eventually form in the waveform. 

The result is illustrated in Fig. 1-1, which is a sequence of time waveforms 
(measured on the left, computed on the right) of a very intense noise pulse as 
it propagates down a plane wave tube.^’^ In the first stage of the distortion 
process, many shocks form as a result of the overtaking of troughs by peaks; see 
the 12 ft waveform in Fig. 1-1. After most of the shocks have formed, 
the distortion process continues but in a qualitatively different way. 

Because each shock has its own particular propagation speed, which is 
determined by its peak and trough pressures, the shock positions in the 
waveform slowly change as the wave continues to propagate. As the shocks 
advance or retreat relative to their neighbors, merging takes place: the 

big shocks tend to "eat up" the little shocks. The noise waveform becomes 
simpler in appearance because the number of zero crossings is reduced. 
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FIGURE 1-1 

NOISE PULSE 1 AT VARIOUS DISTANCES 


(Taken from Ref. 6; 1 ft = 0.3048 m) 









For example, in Fig. 1-1 the propagation from 12 ft to 84 ft is marked by a 
30% reduction in number of zero crossings. In summary, the distortion 
of the noise waveform has two relatively distinct stages. In Stage I 
shocks form but the number of zero crossings does not change very much. 

In Stage II shock merging occurs; accordingly the number of zero crossings 
is reduced. 

The two stages also have, markedly different properties in the frequency 

domain. The formation of shocks (Stage I) is accompanied by a strong 

growth of the high frequency end of the spectrum. In Stage II the 

merging of shocks, which reduces the average time between zero crossings, 

is associated with appreciable growth of the low-frequency end of the spectrum. 

(There is, of course, some low frequency energy produced during Stage I as 

a result of intermodulation distortion. The difference frequency component 

is normally generated much less efficiently, however, than the sum-frequency 

and harmonic components. Distortion in Stage I is therefore manifested 

primarily by high frequency development.) The growth at both ends of the 

spectrum comes, of course, at the expense of the spectral middle. Figure 

1-2, which shows the results of an experimental measurement made by 

6 7 

Pestorius with cw noise, * illustrates the phenomenon. If either of 
the downstream spectra (48 ft or 72 ft) is subtracted from the source 
spectrum (0 ft), the result is the apparent attenuation as a function of 
frequency. The attenuation is negative at low and high frequencies. 

In the midfrequency range the attenuation is much greater than it would 
have been had the noise been small-signal. 

Note that Stage II in the distortion process is not realized unless 
the noise has a very high intensity. At lower intensities, as in the 
Pernet-Payne experiments, ’’steepening" type distortion (Stage I) develops 
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FREQUENCY - Hz 


FIGURE 1-2 

MEASURED NOISE SPECTRUM AT VARIOUS DISTANCES 

(Taken from Ref. 6; 1 ft * 0.3048 m) 
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relatively slowly. Ordinary absorption then has ample time to exercise 

its smoothing effect on the waveform. When this happens, the distortion 
process slows down and eventually stops while still in Stage I. Stage II 

is never reached. 

What modifications to our picture of noise distortion might the 

outdoor environment be expected to introduce? The most noticeable change 

is brought about by the shift from plane to spherical waves. The amplitude 

reduction caused by spherical spreading slows down the distortion process 
13 

considerably. Thus Stage I, which lasted only a few feet in Pestorius’s 
plane wave experiments, is expected to stretch out over a much greater 
distance when the noise is a spherical wave. As we have already seen, 
ordinary absorption opposes the steepening distortion process directly. 
Absorption and spherical spreading thus work in combination against the 
realization of State II distortion. In fact, except for exceedingly 
powerful noise sources, one may expect spectral distortion of outdoor 
noise to be limited to high-frequency buildup. 

Another feature of the outdoor environment that might be expected to 
affect distortion of finite-amplitude noise is random inhomogeneity of the 
medium (random variation in wind and temperature). Our experience in Phase I, 
however, was that while random inhomogeneity definitely affected the 
instantaneous time waveform of the received signal, it had little or no 
residual effect on the long term average spectral content of the signal. 
Random inhomogeneity would not be expected to affect finite-amplitude 
noise any differently. 

The attenuation properties of nonlinear ly generated high-frequency 
noise may be expected to differ from those of small-signal noise of the 
same frequency. In studies of plane periodic waves of finite amplitude in 
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dissipative media (see, for example Refs. 3, 5), it has been found for a wide 
variety of cases that at great distance the second harmonic component 
decays asymptotically as exp ( -2a^x) , not as expC-a^x), where and 02 
are the small-signal attenuation coefficients at the fundamental and second 
harmonic frequencies, respectively. Since generally 2a^ < c^t this means that 
the nonlinearly generated second harmonic component decays more slowly 
than might have been expected from linear theory . Moreover, no distance 
is ever reached at which the second harmonic does decay as expC-c^x), 
i.e., as a small signal. This peculiar behavior is a result of the way 
the second harmonic sound is generated. It is produced not back at the 
source but by nonlinear interaction that takes place over the entire 
travel path between source and observation point. Furthermore, the 
third harmonic component is found to decay asymtotically as exp (-3a^x) , 
not exp ( -a^x) , and the higher harmonics behave similarly. In other words 
the higher harmonic sound never establishes a propagation region in which 
traditional small-signal absorption laws take over. 

How does this peculiar propagation behavior apply to aircraft noise? 
Although the decay laws for spherical waves are different from those for 
plane waves (see Appendix A and Eqs. 2-3 and 2-4 in Chapter 2), it still 
turns out that nonlinearly generated higher harmonics decay less rapidly 
than small-signal sound of the same frequency. This means that although 
high frequency aircraft noise generated at the source decays according to small- 
signal laws, high frequency noise generated nonlinearly in the course of 
propagation decays less rapidly. As distance from the source increases, 
at first the small-signal portion of the high frequency noise masks the 
nonlinearly generated portion. At sufficiently great distances, however. 


11 


only the nonlinearly generated portion survives, and it does not decay at 
the rate expected from small-signal theory. Whether this effect is important 
in practice depends on the distance and amplitude at which the nonlinearly 
generated portion becomes dominant. It does seem clear, however, that the 
commonly held belief that nonlinear effects in aircraft noise may be 
Ignored past a given distance is at best an oversimplification. 

Although this report is primarily devoted to experiments on high 
intensity noise, some data on intense tones are also included. As already 
mentioned, Phase I, which was devoted entirely to tone measurements, is 
reviewed in Chapter 2. Theoretical analyses to support this review 
are given in Appendix A. Some tone measurements were also made during the 
course of Phase II. Data in Appendix C show how SPL of a tone fluctuates 
with time. Finally, one of the sources used for the noise experiments, 
the Applied Electro Mechanics (AEM) 20 Driver Loudspeaker Array, was, at the 
time of our experiments, being evaluated by the Electrical Engineering 
Department for possible use by the Coast Guard as a fog signaling device. 

We carried out tests with high- intensity tones as a contribution to the 
evaluation process. Since the measurements are closely related to those 
made with noise, the results are included in Appendix D. 

The theoretical models used to obtain predictions for comparison 
with the experimental data are described briefly in Chapter 3. Chapter 4 
is devoted to apparatus and procedure. Eleven noise propagation experiments 
were completed. The data, along with linear theory predictions for comparison, 
are presented in Chapter 5. The results are analyzed, interpreted, and discussed 
in Chapter 6. In Chapter 7 an application of the measurements to jet noise is 
given. Here the most important general technological conclusion of the project 
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is drawn: nonlinear propagation distortion is apparently important for 

the noise produced by many jet aircraft. Chapter 8 is a summary. 

Certain ancillary results deemed important enough to report are 
presented in the appendices. As already mentioned, the material in 
Appendices A, C, and D pertains to tones. Appendix B gives a theoretical 
analysis of the amplitude density of finite— amplitude noise. It is 
shown that, prior to shock formation, the amplitude density does not change 
while the wave undergoes nonlinear propagation distortion. 
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CHAPTER 2 


REVIEW OF PHASE I 

In Phase I the outdoor propagation of high- intensity tones was 

studied. The purpose was to determine the extent to which the outdoor 

environment, especially random inhomogeneity of the medium, affects nonlinear 

propagation distortion. Most previous experiments on finite-amplitude 

propagation in air had been carried out indoors under highly controlled 

laboratory conditions, either in tubes (see, for example. Ref. 3, 5, 

and 6 and the earlier works referred to therein) or in anechoic 
X6~18 

enclosures. The outdoor tone measurements were also intended to 

serve as a stepping stone to similar experiments with intense noise 
(Phase II). 


A. EXPERIMENT DESIGN, APPARATUS, AND THEORY 
1. Experiment Design 

Propagation distortion in a spherical wave produced by a monochromatic 

source of radius r Q depends on the source pressure amplitude P^q» source 

frequency f »= to./ 2 tt, propagation distance t, coefficient of nonlinearity 3 

and ordinary absorption of the medium. When ordinary absorption is not 

important, the distortion depends only on the following dimensionless 

19 

combination of factors: 


d - Sekr ln(r/r ) , 

o o 


(2-1) 


which is called the distortion range variable. Here e = p.-Vp c 2 is the 

±u o o 

dimensionless source amplitude, is the static density, c q is the small- 
signal sound speed, and k = go/c q is the wave number. The fact that a is 
proportional to the product of source amplitude and frequency gives us a 
very useful scaling law: if frequency and amplitude are scaled inversely with 
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respect to each other, all other things being held constant, the amount of 
nonlinear distortion In the wave (e.g. , the percent second-harmonic distor- 
tion, percent third-harmonic distortion, etc.) remains the same. When 
ordinary absorption in the medium Is Important, the amount of distortion also 
depends on the dimensionless parameter ar Q » where a is the amplitude attenuation 
coefficient at the source frequency. Because a depends on frequency, the 
.amplitude-frequency scaling law given above is actually only an approximation. 
Even so, the law is still very useful for rough design purposes. 

For more accurate assessment of the ability of a source to produce finite- 

amplitude sound, we developed a graphical method based on the values of the two 

20 

dimensionless parameters f3ekr Q and ar o . See Fig. 2-1, which we call a 
source-frequency level {SFL) chart. The abscissa is ar Q , and the ordinate is 
SFL, whose definition is 

SFL - SPL lm + 20 lo 8io f kHz ^ dB re 20 ^ lPa ~ kHz ” m ^ • ( 2 “2) 

In this formula SPL- is the source level at 1 m (farfield SPL extrapolated 

lm 

to 1 m) arid f ^ is the frequency in kHz. Source-frequency level is a 

21 

convenient measure of 3ekr Q . The lower curve on the SFL chart represents 
source operating conditions that lead to a wave whose calculated shock 
formation distance r (found by setting a = 1 in Eq. 2-1) Is the same as the 
distance r beyond which the fundamental is predicted to decay as a small 
signal. The upper curve is for a much stronger wave, one for which the 
calculated well-formed sawtooth distance & (found by setting o » 3 in Eq. 2-1) 
equals r max * (Formulas for r, / r", and r max are given, for example, in Ref. 19.) 
The SFL chart is used as follows : A source whose operating point falls below 

the lower curve is classified as weak in the sense that its radiation Is 
subject to only weak nonlinear effects. Shock formation Is not expected 
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SFL - dB re 20/z Pa - kHz 







nor much, if any, extra attenuation of the fundamental. Harmonic 
distortion may. develop, the amount depending on the proximity of 
the operating point to the lower curve. Conversely, if the operating point 
lies above the upper curve, the source is classified as strong. Rapid 
distortion, shock formation, the development, of a full-fledged sawtooth, 
waveform, and significant extra attenuation of the fundamental is expected. 
Finally, an operating point lying between the two curves denotes a 
moderate source because it has only moderate ability to generate finite- 
amplitude sound. The various lettered points plotted on the chart are 
operating points of experiments carried out during Phase I. The intensity 
of nonlinear effects observed in these experiments correlated very well 
with the positions of the points relative to the two curves. The SFL 
chart therefore proved to be a valid design tool.* 

2. Apparatus and Measurement 

The experimental facility is now described briefly. To avoid inter- 
ference effects due to ground reflection, we used a vertical propagation 
path that was parallel to a radio tower whose height is 85 m. The source was 
on the ground. After passing through a relatively short nearfield, the waves 
spread spherically over the remainder of the propagation path. The receiving 
microphone was mounted at the end of a 2.6 m boom, which was carried by 
the tower elevator. See Fig. 2-2a. A detailed block diagram of the experiment 
is given in Ref. 1. Alternatively, because the transmit and receive systems 
were similar to the ones used in Phase II, the reader may simply refer to 
Fig. 4-1. 


*To apply the SFL chart to fresh water at 20°C, use 1 pPa-kHz-m as the 
reference for SFL and increase the numbers on the ordinate scale in Fig. 2-1 
by 113.2 dB, e.g., 150 dB becomes 263.2 dB. To allow for salinity and 
temperature changes, make use of information given in Table I of Ref. 21. 
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Three different sources were employed. A 7— element array of * 
exponential horns (Fig. 2— 3a) was used for weak— wave and moderate-wave 
experiments; the horns were driven by JBL 375 horn drivers. The array is 
shown on its mount, a heavy indexing table, in Fig. 2-2b (the baffle visible 
in the picture was not used in any of the propagation experiments) . 

Slightly more intense waves were radiated by a 10-horn array (Fig. 2-3b); 
see, for example, point C in Fig. 2-1. To produce truly strong waves, 
however, we had to construct a siren (Fig. 2 -A) .' With it, experiments with 
operating points in the strong wave region (Fig. 2-1) were possible. 

It was found early that random inhomogeneity of the medium, presumably 
due to velocity and temperature variations, caused fluctuations in the 
received signal. Not surprisingly, the fluctuations increased with 
propagation distance. At the same time, however, it was found that long 
term averaging (the averaging time of the analyzer was selectable) 
provided SPL data that agreed well with predictions based on homogeneous 
media theory. The required averaging time varied from 0.1 sec close to 
the ground to a maximum of 100 sec at the greatest distances. An important 
measurement that could not be averaged was waveform. Waveforms jwere 
recorded with the oscilloscope camera, but they are only samples at 
specific instants. A discussion of fluctuations in Phase II is given in 
Chapter A of this report, and some pure tone data showing the increase in 
fluctuation with receiver height are presented in Appendix C. 

3. Theory 

Three different methods were used to obtain theoretical curves to 

which the experimental data could be compared. First, for weak-wave 

experiments, a perturbation solution of Burgers’ equation for spherical 

2 

waves was obtained. All terms through fourth order were found and also 
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some of the fifth-order terms. Although Burgers T equation applies only to 

2 

a medium whose attenuation coefficient a is proportional of f , the 
atmosphere has this property over limited frequency ranges. The perturbation 

solution is given in Appendix A. Second, an amplitude decay equation that 
had been postulated several years ago for the fundamental component in a 

2 

sawtooth wave was solved (again, see Appendix A). Spherical spreading, 

atmospheric attenuation, and losses due to nonlinear effects (pumping 

of energy to the higher harmonics) are accounted for in this equation. 

Finally, Pestorius's computer algorithm,^ generalized to apply to spherical 

waves in the atmosphere, was used to generate theoretical predictions for 

several of the experiments.^ The algorithm is basically the same as the 

one discussed in Chapter 3. Whenever predictions were computed for a 

given experiment, all the necessary atmospheric attenuation coefficients 

24 

were calculated by the method described in the new proposed standard. 

From the perturbation solution, an important conclusion may be drawn 
about how the harmonic distortion components decay with distance. The 
leading term in the expression for p£, the amplitude of the second harmonic 


component, is (see Eq. A-14 or Refs. 5 and 27) 


fckr o r o -4a (r - r ) _ , ^ 

p 2 = p 10 ~2~~ e ° X 22 (o V “ r) 


(2-3) 


where I ^ is the integral 


’22 


/ 


r 2a(r' - r ) 
e o 


dr’ 


34 

which can be expressed as the difference of two Ei functions. Note that 
1^2 increases monotonically with distance. Although the second harmonic 
component has some resemblance to a small-signal wave of frequency 2f 


* D. T. Blackstock, unpublished research (1971) . 
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(notice the factors denoting ordinary spherical spreading and exponential 
attenuation in Eq. 2-3), it really behaves differently. After first 
increasing with distance, p^ reaches a maximum and then decreases. The 
decay region is particularly interesting. Because of the monotonically 
increasing property of I ,,, t * ie second harmonic never reaches a point at 
which its decay approaches that of a small signal. It always decays less 
rapidly than a small-signal wave of the same frequency. In particular, 
at great distances, the formula for p£ is 

Bek / r o\ -2a(r - r ) /o 

pa ~ p io vry e ' ° • (2_4) 

which deviates markedly from what would be expected on the basis of linear 

2 

theory. It will be seen that P 2 is proportional to p^ , where 

r, m D ( r /r)e"* a ^ r ” r o^ is the local amplitude of the fundamental. 

*1 *i0 o 

Similarly it can be shown that the third harmonic amplitude never 

behaves as a small signal of frequency 3f; in fact, at great distances p^ 

3 

is proportional to . The deviation from small-signal behavior therefore 
increases with harmonic number. The application of this result to air- 
craft noise was discussed in Chapter 1. 

Although the theoretical models proved generally successful in 
explaining the experimental measurements, there was one recurring problem: 
how to account for nearfield distortion. Our theoretical models are for 
true spherical waves, which emanate from a source of radius r Q . Although 
sources of the type used in our experiments do produce spherical radiation, 
there is first a nearfield. Spherical Bpreading begins at approximately 
the Rayleigh distance (& q » S/A, where S is the radiating area of the 
source and A is the wavelength). Because some distortion takes place in 
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the near field, the wave is not monochromatic, as assumed in our simple 
theory, when spreading begins. A rigorous analytical treatment of finite- 
amplitude distortion in the nearfield was out of the question. Our near- 
fields were too complicated. At the same time, nearfield distortion was too 
large to be ignored. We therefore accounted for it empirically. One ,empiri- 

19 

cal method that has been used successfully in the past is to replace the real 

source and its niearfield by an equivalent spherical source that matches 

directly to an inward extension of the farfield. The radius r of the 

o 

equivalent source is taken to be a fraction of the Rayleigh distance. 

The fraction is to be chosen so that the equivalent spherical wave distorts 

as much in traveling from r to R as the real wave does in traveling 

o o 

through the nearfield. In the past, giving the fraction a value in the 

19 

range 1/3 to 3/4 has worked well. An arbitrary choice, such as 

r Q * R q / 2, is useful for design purposes, for example, for specifying an 

operating point on the SFL chart. For application to a particular experiment, 

however, once the data have been taken, the values of p^ and p^ measured 

at a single point In the farfield can be used to compute the value 

of for that experiment. We used this procedure to obtain values of r^ 

for our theoretical predictions. 

The equivalent spherical source model has several drawbacks, however. 

The near field phase changes, which can be large, are ignored in the model. 

The fact that r Q is frequency dependent makes the model difficult to apply 
when the source signal itself is distorted or broadband.* Finally, in 
practice we found that r Q varies with source amplitude as well as frequency. 

*In fact, the model proved to be of little use in Phase II, where the source 
signal was noise. 
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Therefore, despite its conceptual simplicity, the model leaves much to be 
desired. It is useful for design calculations and for providing an initial 
match between theory and measurement for a specific experiment. It is not 
however, worthy of blind faith. 

B. EXPERIMENTS 

1. Weak Waves 

Figure 2-5 shows the results of an experiment done with the 7- 

horn array at 8.25 kHz. All data points shown are in the far field; the 

Rayleigh distance was 3.2 m. On the SFL chart (Fig. 2-1) the operating 

point for the experiment is A' (SFL = 160 dB, ar^ = 0.012, where for this 

purpose the value R / 2 is used for r ) . Although A 1 falls slightly above 

o o 

the lower curve, for simplicity we call this a weak-wave experiment. 

For theoretical calculations the equivalent spherical source model was not 
used to compensate for nearfield distortion. Instead, r^ was taken to 
be the initial measurement distance (6.1 m) . The distorted signal measured 
at this point was specified as the boundary condition for the perturbation 
solution of Burgers’ equation.* The three curves pertaining to p ^ are 
explained as follows: The dotted curve represents the component of the 

second harmonic generated nonlinearly in the region r > r Q (Eq. A-9a) . 

The solid curve represents the small-signal decay of the second harmonic 
component Initially present at i = r^. When the two components are com- 
bined according to Eq. A-21 (in this case the phase angle <f> was approximately 
0°), the result is the dashed curve, which accounts very well for the 
measured data. As for the fundamental, the solid curve represents the 
linear theory prediction. The dashed curve includes the reduction of the 
fundamental caused by its nonlinear interaction with the Initial second 

* Burgers’ equation is suitable for predicting p 1 and p^ because, for the 
conditions of the experiment, = 4a^. ^ 



RECEIVED SPL - dB re 20,4 Pa 



FIGURE 2-5 

PROPAGATION DATA, WEAK-WAVE EXPERIMENT 
SPL = 141.5 dB, f = 8.25 kHz, TEMP = 27*C 
RH = 72%, a = 0.0076 Np/m. 

(POINT A' IN Fig. 2-1) 
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harmonic (see Eq. A-20) . The dashed curve seems to fit the data slightly 
better than' the solid curve. In summary, second-order perturbation theory 
provides an excellent explanation -of both propagation curves. Moreover, 

the mildness of the nonlinear effects observed in this experiment correlates 
well with the relative position of point A’ on the SFL chart. 

Notice how the second harmonic signal generated nonlinearly in the 

region r > r Q (dotted curve) eventually overtakes and masks the second harmonic 

that was Initially present at r . Although the latter was itself non- 

o 

linearly generated In the region r < r^, it is small enough in amplitude 
to decay as a small signal in the region r > (note that the decay factor 

is e a ^ r r o\ not e ar ) . Considering as the effective source radius, 

therefore, we see that the behavior of the data in this experiment 
supports an assertion made in Chapter 1: When a nonlinearly generated 

second-harmonic sound is masked by a small signal of the same frequency, 
the propagation will eventually be dominated by the nonlinearly generated 
component because of its less rapid decay. Unfortunately, our range of 
measurement distance did not extend far enough to test the asymptotic 
formula Eq. 2-4. 

2 . Moderate Waves 

The results shown in Figs. 2-6 and 2-7 are for an experiment done with 
the 10-horn array. The Rayleigh distance was 3.6 m, the frequency was 
6.6 kHz, and the source level was 146.5 dB. The experiment is represented 
on the SFL chart (Fig. 2-1) by point C (SFL = 163 dB, ar Q * 0.0096), which is 
right in the middle of the moderate-wave region. The theoretical predictions 
were obtained by using the computer algorithm. The measured waveform at 
6.1 m (Fig. 2-7a) was used as the initial signal for the computation. 
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PROPAGATION DISTANCE - m 


FIGURE 2-6 

PROPAGATION DATA, MODERATE-WAVE EXPERIMENT 

(POINT C ON Fig. 2-1) 
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EXPERIMENTAL 

WAVEFORM 


COMPUTED 

WAVEFORM 


SPECTRA 






X: EXPERIMENTAL 
BAR: COMPUTED 
f = 6.60 kHr 


FIGURE 2-7a 

WAVEFORMS AND SPECTRA, MODERATE-WAVE EXPERIMENT 
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EXPERIMENTAL 

WAVEFORM 


COMPUTED 

WAVEFORM SPECTRA 





X: EXPERIMENTAL 
BAR: COMPUTED 
f = 6.60 kHz 


FIGURE 2-7 b 

WAVEFORMS AND SPECTRA, MODERATE-WAVE EXPERIMENT 







During the course of the measurements the atmosphere was very calm. 

There was no wind at all on the ground and the anemometer read only 0-4 km/hr 
at the top of the tower. 

Nonlinaar distortion in this experiment is seen to be appreciably 
greater than in the previous one. The fundamental clearly suffers extra 

attenuation, albeit not large. Shock formation is indicated by the 
crowding together of the propagation curves in Fig. 2-6 and is confirmed 
by the waveforms in Fig. 2-7. In fact, the waveforms at 12 m and 21 m 
have a nearly sawtooth shape. When this experiment is contrasted with the 
previous one, it is interesting to note that the large increase in nonlinear 
effects is due to an increase of only 3 dB in SFL (compare the positions 
of points A T and C in Fig. 2-1). The relatively sharp onset of finite- 
amplitude behavior was noticed repeatedly throughout the project. 
Nonlinearity in spherical waves seems to be almost a threshold phenomenon. 

An increase of 3 dB in source level can have either a small effect or a 
large effect, depending on the starting point. 

Except at the greatest distances, the waveforms In Fig. 2-7 are 
asymmetric. The peaks are sharp and the troughs are round. Deviation of 
this sort from the ideal sawtooth shape is probably caused by diffraction 
effects, which are not included in the simple spherical wave models dis- 
cussed thus far. There are at least two sources of diffraction in our 
experiments. First there is piston type diffraction. If the mouth of a 
horn Is likened to a radiating piston, the following result from linear 
piston radiation theory may be used: on axis the farfield signal Is 

proportional to the time derivative of the piston velocity.* In this 

* In Chapter 3 we discuss the frequency domain version of this result, 
which is that the farfield spectrum has a 6 dB/octave boost relative 
to the mouth spectrum. 
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experiment the mouth signal is not a sinusoid. In traveling through the 
horn, the wave suffers some finite-amplitude distortion. It arrives at 
the mouth with a waveform similar to the last measured waveform in Fig. 2-7b. 
Given this shape, differentiation will produce a waveform that is asymmetric 
like the one at 6.1 m in Fig. 2-la, This effect was observed and reported 
20 2 

early in Phase I. * Another diffraction effect, which also produces 

the same kind of waveform asymmetry, is associated with narrow beams. 

23 

This effect has been analyzed by the Russian investigators. Although 
our beam does not have quite the same characteristics as the beam they have 
analyzed, our beam was relatively narrow (3 dB beamwidth of about 6° in this 
experiment). Some of the asymmetry evident in Fig. 2-7 may therefore be due 
to narrow beam diffraction. 

Two further aspects of the asymmetry may be noted. First, although 
the asymmetry eventually disappears as propagation distance increases, the 
disappearance is faster in the computed waveforms than in the measured 
waveforms. The computed waveforms become symmetric more rapidly because 
(1) no diffraction effect is included in the algorithm to perpetuate the 
initial asymmetry, and ( 2 ) ordinary nonlinear steepening tends to wash 
out the asymmetry. The fact that in the measured waveforms the asymmetry 
seems to first increase (or at least hold its own) with distance supports 
the speculation that narrow beam diffraction may have been important 
while the beam diameter was still small. Second, as noted in the 
footnote on the preceding page, the differentiation operation, which accounts 
for the asypnetry (no matter whether the asymmetry is due to source diffraction 
or beam diffraction), gives the higher harmonics a boost. This may help 
explain why the second harmonic data lie somewhat above the computed 
propagation curve in Fig. 2-6. In fact, at 12 m the separation between 
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data points for the fundamental and second harmonic signals is slightly 
less than 6 dB. A separation of less than 6 dB is not expected for an 
ordinary sawtooth wave, but as Fig. 2-7 shows, the waveform at 12 m is a 
peaked sawtooth. 

Finally, the discrepancy in Fig. 2-7 between predicted and measured 
waveforms at 76 m, and to a lesser extent at 61 m, is probably due to 

inhomogeneity of the medium. Although each data point in Fig. 2-6 represents a 
time average measurement, the waveforms in Fig. 2-7 are single samples. 

As the fluctuations increase with propagation distance, so does the possibility 
of recording a sample that varies considerably from the average. 

3. Strong Waves. 

Propagation data for one of the siren experiments are shown in Fig. 2-8. 

The short conical horn (diameter 0.152 m) fitted to the end of the siren 

(see Fig. 2-4) was designed to convert the siren from a ring source into 

a piston-like source. The operating conditions for the experiment described 

here, f=6.33 kHz, SPL^ m = 150 dB, and a = 0.0070 Np/m — define point G on 

the SFL chart (Fig. 2-1). The transmitted sound had a wide beam, approximately 

40° between 3 dB down points, and a short nearfield (the computed value of 

R was 0.33 m) . 
o 

The sound generated in this experiment was very intense. The . 
magnitude of the intensity may be judged in several ways. First, the 
crowding together of the propagation data for the fundamental, second 
harmonic, and third harmonic, even at distances less than 1 m, suggests 
that the wave formed a sawtooth very close to the source. Indeed this 
was the case, as oscillograms of the waveform show.'*’ Second, the propagation 
data for the fundamental deviate markedly from the linear theory curve . 

At 30 m the deviation, or extra attenuation, is 10 to 11 dB. Third, from 
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about 3 m on, the data are only 1 to 2 dB less than their predicted saturation 

values. At saturation the received level is at an absolute maximum value 

19 24 

that cannot be exceeded no matter how high the source level is raised. * 
Therefore, raising the source level of the siren from 150 dB to, say, 156 dB 
(a quadrupling of the source power), would not be expected to change the 
propagation data much from that shown in Fig. 2-8. 

Although not shown in Fig. 2-8, theoretical predictions obtained by 
means of the computer algorithm provided a reasonably good fit to the 
propagation data. Predictions based on a simple amplitude decay rate model, 
explained in Section II of Appendix A (where a second siren experiment is 
also described) , were found to agree very well with the data for the 
fundamental. See Eq. A-28, from which it is a simple matter to obtain 
predictions of both extra attenuation and saturation amplitude. 

The sawtooth waveforms in this experiment were also observed to be 
asymmetric.^ The asymmetry was strong up close to the source but tended 
to die out monotonically and more rapidly than in the 10-horn array 
experiment (Fig. 2-7). This behavior is consistent with the difference 
between the sources in the two experiments. Because the siren produced 
a very broad beam, asymmetry was limited to that caused by source diffraction. 
Because of the greater intensity in the siren experiment, the asymmetry 
was more quickly overcome by nonlinear steepening. 

C. Summary 

In summary, the tone experiments covered a wide range of wave strengths, 
from very weak to very strong. The classical nonlinear propagation phenomena 
that are well known from indoor experiments — harmonic distortion, shock 
formation, and extra attenuation (to the brink of saturation) — were observed 
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FIGURE 2-8 

PROPAGATION DATA, STRONG-WAVE EXPERIMENT 
(POINT G ON Fig. 2-1) 



in the outdoor environment. Moreover, provided sufficient averaging time 
was allowed for the measurements, theory and measurements were found to be 
in good agreement even though the theory was for a homogeneous medium. 
Asymmetric waveforms qualitatively attributed to diffraction effects were 
ob served. 
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CHAPTER 3 


THEORY 


In this chapter the theory of freefield propagation of finite- 
amplitude noise in the atmosphere is considered. Although the atmosphere is 
actually inhomogeneous, our experience in Phase I showed that, for our 
experimental conditions, the assumption that the medium is quiet and 
homogeneous works well if sufficient signal averaging is employed. In our 
discussion we therefore consider freefield propagation in a homogeneous, 
quiet medium. 

A. LINEAR THEORY 


Let F(r ,w) be the noise power spectrum at range r > R , where the 
o o — o 

Rayleigh distance R q marks the beginning of the transmitter farfield. 

If the subsequent propagation is dominated by spherical spreading and 
atmospheric absorption, the power spectrum at range r >_ r Q is given by 

r 2 e [- 2 a(m)(r - r Q )] 

F(r,w) = — £2 F(r o ,w) , (3-1) 

where a (to) is the atmospheric absorption coefficient at angular frequency to. 

In this report absorption coefficients are calculated according to 

22 

the new proposed standard. We shall frequently use the term "linear 
theory" to characterize computations based on Eq. 3-1. 

Nothing has been said so far about the source of the farfield 
spectrum F(R q , to). We take the source to be a circular horn or aperture. 

The geometry is shown in Fig. 3-1. The wavefront is assumed 
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FIGURE 3-1 

GEOMETRY OF THE DIFFRACTION OF SOUND 
BY A CIRCULAR APERTURE 
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to be plane in the aperture; the power spectrum there is given by 
F(OiU>). The spectrum at any point r on the axis in front of the aperture 
is 


F(r ,«) - | e -<J fc+o)r - e _(jk+a)r l | 2 F<0, 


0 )) » 


(3-2) 


1/2 


where J“(-l) and r^-i/r z + a z is the distance from the aperture 

edge to the observation point (see Fig. 3-1). In the farfield (r»a) » , 

r^ ■ r + a 2 /2r and Eq. 3-2 reduces to 

F(r,u) - ) (l + |z) e‘ 2or F(0,«O , 


/R (w)\* . 

■- (-V - ) e F<0 >“ ) 


(3-3) 


since a/k «1 for f <1 MHz. In terms of the aperture area S=7ra 2 , 
the Rayleigh or farfield distance is 


V“> 


So) 

2irc 


o 

Because the farfield distance increases with frequency, the farfield 
spectrum has a high frequency boost compared to the source spectrum. 
Provided the observation point is not deep in the farfield, F(r ,a)) oc wF(0, 
i. e., the boost is 6 dB/octave. The farfield distance for the entire 
spectrum is R o (io m ) where corresponds to the highest frequency of 
interest in the source spectrum F(0 r w). For other sources of interest, 
e.g., a rectangular aperture, Eqs. 3-3 and 3-4 may be used by taking S 
to be the radiating area of the source. If the source is an array of 
apertures, S is the total area of all apertures. We shall refer to com- 
putations based on Eqs. 3-2 and 3-3 as "linear diffraction theory." 


u) 


(3-4) 


> 
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B. NONLINEAR THEORY 


In the previous section we presented results valid in both the trans- 
mitter nearfield and farfield. As noted in Chapter 2, there is no well- 
developed theory for the nearfield of finite-amplitude waves. In analyzing 
propagation of finite-amplitude waves, therefore, we have to restrict our 
attention to the farfield problem. 

Our approach to the problem of finite-amplitude spherical noise 
waves is similar to that of Pestorius, who developed a computer algorithm 
to treat plane waves of noise. ^ The analysis we present here is valid 
prior to shock formation, although it can easily be extended beyond by using 
the relations of weak-shock theory.^ The plane wave algorithm is based on 

the solution of the nonlinear wave equation for plane progressive waves in 

2 

a lossless gas. The equation is c u + c u - $uu = 0, where x is 

o x o t t * 

propagation distance and t Is time, or, if the retarded time t ,=, t-x/c o 
Is introduced, 

c 2 u - 3uu , = 0 . (3-5) 

ox t 

This equation is the mathematical embodiment of the distortion process, 

wherein the compression phases of a wave tend to steepen; see 

Chapter 1. For a boundary condition of the form 

u | = g(t) = g(t') , (3-6) 

! x =o 


the exact solution of Eq. 3-5, called the Earnshaw solution, is 


where 


t' 


U = g(T) , 
Sxu(t) 


T - 


(3-7a) 

(3-7b) 
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The parameter t is the time base for the source signal, that is, 

T“t'"t when x“0. Although it is difficult to obtain analytical 
results from the Earnshaw solution, the solution may easily be 
implemented as a computer algorithm. The waveform at any desired 
distance x is computed simply by distorting the time base according to 
Eq. 3-7b. However, ordinary absorption, which is frequently important, is 
not included in the Earnshaw solution. In order to account for ordinary 
absorption, we assume that even when the medium is lossy, Eqs. 3-7 may be 
used to propagate the wave a suitably small distance Ax. After a transforma- 
tion to the frequency domain, absorption corrections appropriate for the 
distance Ax are applied. There follows a transformation back to the time 
domain, where another short propagation step is made. The process is 
repeated until the desired propagation distance is reached. 

For spherical waves, the computation is much the same. The wave 
equation for spherical waves reduces to Eq. 3-5 if one makes the changes 
of variable 


r 

w =» — u , 

r 

o 

and 

x = r ln(r/r ) 
o o 

In terms of distortion, spherical wave propagation from r to r + Ar is 


equivalent to plane wave propagation from x to x + Ax if Ax is given by 


(3-8a) 

(3-8b) 


Ax - r ln[l + — ] 
o r 


(3-9) 


A block diagram of the computational procedure for spherical waves, which Is 
based on Eqs. 3-7 and 3-8, is shown in Fig. 3-2. It is assumed here that we 
wish to propagate the wave from r = r^ to r =* r^, where r^ >_ r 0 * Ttie input 
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FIGURE 3-2 

BLOCK DIAGRAM OF THE COMPUTATIONAL PROCEDURE 
FOR SPHERICAL WAVES 

(NOTE THAT IF SHOCK FORMATION IS POSSIBLE, THE 
k-2 BLOCK LABELED PROPAGATE COULD BE REPLACED 

WITH THE WEAK SHOCK PROPAGATION MODEL OF Ref. 6) 













* 

waveform w is specified at r ■ r , as is the incremental step size Ax. 

A check is performed to determine whether propagation over Ax will exceed the 

desired range r^. If so. Ax is set to a value such that r + Ar ■ r^. If not, 

the current value of Ax is used. The equivalent spherical wave step size Ar 

is computed, and the wave is propagated the equivalent plane wave distance 

Ax according to Eqs. 3-7. We next transform to the frequency domain in order 

to make absorption corrections. After a transformation back to the time 

domain, the current value of r is incremented. When the desired propagation 

distance is reached (r = r^) , the correction for spherical spreading is 

made (see Eq. 3-8a) and the computation is complete. 

In order to use the computational procedure described in the previous 

paragraphs, we must have a representation of the source waveform u(r ,t). 

o 

One approach would be to record a sample waveform of the noise at r ■ r Q 
and use this as the source waveform. There are practical problems associated 
with this approach, however. First, the sample length must be sufficiently 
long to reflect both the statistical and spectral properties of the noise. 

In an outdoor environment, signal fluctuations caused by the randomness 
of the medium may make the required sample length prohibitively long. 

Second, recording equipment is required. Finally, even in field studies 
in which recording equipment is used, the results are almost always 
reported in terms of spectra, not waveforms. Developing an input 
waveform from a spectral measurement therefore seemed to be an attractive 
alternative to recording the waveform directly. 

A method for constructing a suitable source waveform s(t)"u(r Q ,t) given 

*The step size Ax is typically chosen as some fraction (say 0.1) of the 
shock formation distance appropriate for a lossless medium. See Ref. 6 
for further details. 
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just the power spectrum F(r Q ,w) is given in Ref. 25. We start with a 

line spectrum of closely spaced components at the harmonic frequencies nAf. 

The envelope is a replica of the given average power spectrum F(r o> u)), 

scaled so that the mean square values of the two signals are equal (see 

Fig. 3-3). If it is assumed that the given noise is Gaussian, the phase of 

each component in the given noise is a random variable $ whose probability 

25 

density is uniform over the range -ir <_ <f> Tr * In order that u(r Q ,t) have 

the same property, we assign a random phase with a uniform probability 

density to each of the harmonic components. Now having a complete frequency 

domain description of the desired noise, we use an FFT 1 operation to 

obtain the time waveform u(r ,t). The difficulty with this method is, of 

o 

course, that there are an Infinite number of signals u(r Q ,t) with the same 
power spectrum, i.e., there exist an infinite number of ways to assign 
the random phase $>. Clearly the success of our method for representing 
u(r Q ,t) depends upon the sensitivity of the computations to the particular 
function u(r Q ,t) chosen. In a study of the interaction of an intense 
tone with low level noise, Webster and Blacks tock found that the method 
described here worked quite well. For the case of intense noise, however, 
the validity of the model rests on comparison with the experimental data. 
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F(f) FFT’ 1 
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FIGURE 3-3 

CONSTRUCTION OF THE SOURCE WAVEFORM GIVEN THE 

POWER SPECTRUM F(f) 
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CHAPTER 4 


EXPERIMENTAL APPARATUS AND PROCEDURE 


In this chapter we describe the equipment and procedures used to make 
the experimental measurements reported in Chapter 5 and Appendices C and D. 
The experimental arrangement is much the same as that used in Phase I of 
the project. ^ We shall therefore give only a general description of the 
apparatus and procedures and point out the changes Instituted for Phase II. 
For a more detailed account, the reader should consult Ref. 1. 

A. EXPERIMENTAL APPARATUS 

A general schematic of the experimental apparatus is given in Fig. 4-1. 
The electroacoustic source was either a 20-element horn array manufactured 
by Applied Electro-Mechanics (AEM) or a horn array driven by JBL 375 
horn drivers. The AEM array (shown in Fig. 4-1) Is described in detail in 
Appendix D. The JBL arrays were like the 7- and 10-element sources used in 
Phase I (Fig. 2-3), but single-element and 3-element arrays were also used. 
The indexing table on which the source was mounted (see Figs. 2-2b and 2-3) 
could be rotated and tilted for alignment and beam pattern measurements. 

The propagation distance was varied by moving the elevator car, on which the 
microphone boom was mounted, up and down the tower. In addition to the 
traveling microphone (1/4 In.) shown in Fig. 4-1, a monitor microphone 
(1/4 in. or 1/8 in.) was used for measurments near or inside the horn array. 
For all measurements taken with the traveling microphone, however, the 
monitor microphone was removed so that it did not affect the wave field. 

A list of titoe more important equipment used is given below. For a 
more detailed list, see Ref. 1. 

(1) Microphone. B & K type 4136, 1/4 In. The -pressure response is 
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ELEVATOR 



FIGURE 4-1 

BLOCK DIAGRAM OF THE EXPERIMENTAL ARRANGEMENT 













nominally flat from 50 Hz to 50 kHz (1 dB down points) for grazing (90°) 

incidence. If the microphone is used under freefield conditions, the 

manufacturer's rated corrections to be added to the pressure response 

are -0.2 dB for the frequency range 9-14 kHz and +0.5 dB for the range 
26 

19-35 kHz. All other components in the receiving system have a pass 

band at least as wide as 20 Hz to 100 kHz. The microphone is therefore 

the limiting element in the receiving system. 

(2) Microphone. B & K Type 4138, 1/8 in. The 4138 microphone 

response was rated as flat to within ±1 dB from 20 Hz to 100 kHz. No free- 

26 

field corrections were required to 70 kHz. This microphone was used as 
a monitor only. 

(3) Microphone Preamplifier. B & K Type 2619. 

(4) Microphone Preamplifier and Power Supply. B & K Type 2801. 

(5) Microphone Cables. B & K Type A0 0029 (30 m) and AO 0028 (10 m) . 

(6) Pistonphone (microphone calibrator). B & K Type 4220, 

(7) Spectrum Analyzers. B & K Type 2010, H-P Model 3580A and 
G-R Type 1900-A. 

(8) Power Amplifier. CML McCarr 5 kW. 

(9) Current Transformer. Pearson Electronics Model 110. 

(10) Signal Generators. B & K Type 1022 Oscillator or Elgenco 
Gaussian Noise Generator. 

(11) Filters. B & K Type 2112 Audio Spectrometer and K-H Model 
335R or 310-AB Band Pass Filters. 

B. EXPERIMENTAL PROCEDURES 

In this section we describe the procedures used in equipment checkout 
and calibration and also in the actual data acquisition. The procedures 
were the same for both noise and tone signals. 
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1. Equipment Checkout and Calibration 

Prior to and following each of the experiments, the following tasks 
were performed. First, the particular spectrum analyzer used was 
calibrated in both frequency and amplitude. Second, the pistonphone was 
used to make an absolute calibration of the receive system. The calibration 
was done first with the microphone connected directly to the microphone 
power supply and then with the 100 m cable (composed of three 30 m cables 
and one 10 m cable) inserted between the two (see Fig. 4-1) . Comparison 
of the two calibration values indicated whether the cable and its 
connections were sound. If the two calibration values differed by more 
than 0.2 dB (the accuracy of the pistonphone output), cable repair was in 
order. Moreover, comparison of day to day calibration values was made to 
detect possible equipment malfunctions. Again, if deviation of more than 
0.2 dB was observed, each component in the receive system was checked to 
determine the cause. 

2. Data Acquisition 

The procedure was as follows. First, the elevator was moved to the 
desired position. The source was then activated and the received spectrum 
recorded. For the noise measurements, both the H-P 3580A and B & K 2010 
analyzers were used. The H-P analyzer was used to record the continuous 
spectrum, while the B & K analyzer was used to measure the overall sound 
pressure level (OASPL) . As a crosscheck on the measurements, the B & K 
analyzer was also used to record 100 Hz band levels at specific frequencies 
(usually 5, 10, and 20 kHz). These measurements were then compared with 
those recorded by the H-P analyzer. If the two measurements differed by 
more than 1 dB, the measurements were rejected and the cause of the discrepancy 
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was found. The rms input voltage and current to the array were also monitored 
during the measurement period. After all data were recorded at each of 
the desired ranges, the elevator car was moved back to the first measurement 
position. A new spectrum was recorded and checked against the spectrum 
previously recorded there. If the two spectra differed by more than 1 dB, 
the data were rejected and the cause of the difference was identified. 

The system of checks described here proved invaluable in assuring the 
internal consistency of the data. 

Some mention should also be made of the averaging time used for the 
spectral measurements. In addition to the averaging required because of 
the statistical nature of the signal (noise), inhomogeneous medium effects 
(gustlng winds and temperature variations) caused severe fluctuation 
in short-term average levels (see Appendix C) . The averaging time of the 
H-P analyzer could be Increased by increasing the RC time constant of the 
postvideo detector filter (the constant was adjustable via a front panel L 
knob). The ’’maximum smoothing” setting (maximum averaging time) was found 
to be quite adequate when the 100 Hz bandwidth was used. The averaging time 
for the B & K analyzer was adjustable from 0.1 sec to 100 sec. For OASPL 
measurements averaging times of 1 sec to 3 sec were used; for 100 Hz 
band measurements the 100 sec position was chosen. 

C. SUMMARY OF EXPERIMENTAL CONDITIONS 

A catalog of the conditions for the experiments conducted In the 
Phase II research is given in Table 4-1. For the noise measurements the 
spectrum band limits correspond roughly to the most intense portion of the 
spectrum. The source level (OASPL^) is the farfield OASPL extrapolated 
(by spherical spreading) to 1 m. For the tone measurements OASPL^ is 
the extrapolated farfield SPL of the fundamental frequency component. 
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The characteristic ground level wind speeds observed roughly divide 
the experiments into two groups. During experiments JBL 3, 4, and 7 
and AEM 1, 2, 5, 6, 8, and 9 wind speeds generally ranged from 13 to 24 
km/h. For the remainder of the experiments wind speeds were between 7 
and 10 km/h. 

Finally, the experiments were performed at night during the period 


June-September 1977. 



TABLE 4-1 


CATALOG OF EXPERIMENTS 


Experiment Figure 
Number 

OASPL, * 
lm 

dB 

Input Spectrum 
Band Limits, 
kHz 

V £) 

m(kHz) 

Temp, 

(°C) 

. RH 
(%) 

JBL 

1 

5-6 

121 

1-7 

0.53 

(10) 

23 

75 


2 

5-7 

127 

3 - 5.5 

1.06 

(20) 

26 

73 


3 

5-8 

132 

3.5 - 6.6 

3.18 

(20) 

24 

90 


4 

5-9 

144 

3.5 - 7 

7.42 

(20) 

24 

86 


5 

5-5 

140 

2.8 - 5.6 

7.42 

(20) 

26 

73 


6** 









7 

5-10 

142 

5.8 - 7 

7.42 

(20) 

28 

55 


8 

5-11 

141 

3.5 - 6.5 

7.42 

(20) 

27 

65 


9 

5-1 

145 

3-8 

10.60 

(20) 

29 

61 

AEM 

1 

D-5 

127 

5.33 (Tone) 

4.06 

(5.33) 

26 

82 


2 

D-6 , D-7 

143 

5.33 (Tone) 

4.06 

(5.33) 

28 

61 


3 

5-12, 6-2 

132/142 

2-5 

7.60 

(10) 

27 

72 


4 

5-2, 5-3, 5-4 

123/142 

2.8 - 5.6 

7.60 

(10) 

26 

82 


5 

D-l, D-2 

128/146 

1.5 (Tone) 

1.14 

(1.5) 

28 

72 


6 

D-3, D-4 

125/144 

3.55 (Tone) 

2.71 

(3.55) 

31 

57 


7** 









8 

D-9 

127/147 

5.33 (Tone) 

4.06(5.33) 

28 

73 


9 

5-13 

133 

2.8 - 5.6 

7.60 

(10) 

29 

61 


When both low level and high level measurements were done, both source levels 
are given. 

Experiment not completed. 
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CHAPTER 5 


EXPERIMENTAL RESULTS AND COMPARISON WITH LINEAR THEORY 
A. INTRODUCTION 

In this chapter we present the results of the noise measurements. 

Two experiments, JBL9 and AEM4, which exhibit the main features of the 
nonlinear phenomena, are described in detail in the text. The novel 
features of the remaining nine experiments are discussed. In each case 
the measurements are compared with predictions based on linear theory 
(Chapter 3) . A comparison of the data with predictions based on nonlinear 
theory (Chapter 3) is made in Chapter 6. 

Several properties of the spectral information and its presentation 

(see, for example, Fig. 5-1) are common to all the sets of data. First, 

the spectra were obtained with a fixed bandwidth (typically 100 Hz) 

spectrum analyzer (HP-3580A) . As is always the case with such an analyzer, 

measurements within a few bandwidths of zero frequency are not reliable. 

Moreover, wind generated noise often contaminated measurements below 

500 Hz. For these two reasons, the lower limit of useful data Is about 

500 Hz. When a measured spectrum was used as the starting point for a 

theoretical prediction, the lower (suspect) end of the spectrum was 

replaced by a gradual extrapolation to zero. Second, at high frequencies 

the lower (SPL) limit of useful measurements is either 30 dB below the 

middle horizontal grid line (indicated by tick marks) or the absolute level 

SPL - 25 dB (100 Hz band), whichever is higher. The first limit is set by the 

dynamic range of the spectrum analyzer, the second by other instrumentation noise 

in combination with ambient acoustic noise. Third, the distance R (f ) is 

o m 
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meant to represent the nearest point at which all significant spectral 

components radiated by the source have reached the farfield. Since R q 

increases with frequency (see Eq. 3-4), f should be chosen to be the 

m 

highest frequency of interest present in the noise spectrum at the 

aperture, in this case, the face of the transmitting array. In cases * 

where the aperture spectrum had been measured, f was so chosen. If 

the aperture spectrum had not been measured, f was taken to be the 

m 

highest frequency in the spectrum analysis. Fourth, the indicated 

value of OASPL is an effective source level at 1 m. It is equal 

to OASPL | + 20 log 1Q r Q , where r Q is normally the distance (in m) of the 

first measurement point beyond P^f^). We justify the use of spherical 

spreading to extrapolate back to 1 m by assuming that although the noise 

wave distorts as it travels to r , little energy is actually lost. 

o 

That is, few if any shocks form, and ordinary absorption is small over that 
distance. Finally, two different procedures were used in computing the 
predictions based on linear theory, depending on whether the effect of 
diffraction was included. When only spherical spreading and atmospheric 
attenuation were taken into account, the starting spectrum for the calculations 
had to be one measured in the farfield. The spectrum measured at r^ was used. 
This procedure was followed for Figs. 5-1, 5-4 through 5-11, and 5-13. 

When diffraction was also included, as in Figs. 5-3 and 5-12, the starting 
spectrum was one measured inside one of the horns. Ordinary linear horn theory 
was used to extrapolate to the horn mouth (the aperture), and Eq. 3-3 was used 
to bridge the nearfield. Predictions are limited to ranges r r . 

B. EXPERIMENTAL RESULTS 

A typical set of spectral measurements is shown in Fig. 5-1. The 
spectrum of the array input current is shown at the top of Fig. 5-la. 
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Because range was measured from the face of the array, negative values of 
r signify points inside one of the horns. For the acoustic measurements 
the ordinate is the SPL in a 100 Hz band (the indicated bandwidth in all 
measurements reported here is the equivalent filter noise bandwidth) and 
the abscissa is the filter center frequency. The data in the left column 
of Fig. 5-la are for points inside the horn. In the right column 

the measurement point was in the nearfield for most of the spectral components. 
For Fig. 5-lb all the points were in the farfield. In the horn and in 
the nearfield, the growth with distance of the high frequency portion 
of the spectrum is evident. The growth is caused by both nonlinear 
distortion and diffraction. In Fig. 5— lb a comparison Is made of the 
farfield data with linear theory [in this case 9.14 m was used as 
the value of r even though that value is slightly less than R (f )]. 

The linear theory predictions of the most intense part of the spectrum 
(approximately 3 kHz to 7 kHz) agree well with the measurements. At 
higher frequencies however, the data rise above the predicted levels; 
the discrepancy increases with range and frequency. At the 79.3 m point, 
for example, the measured level exceeds the predicted level by 2 dB 
at 10 kHz, 6 dB at 15 kHz, and 15 dB at 20 kHz. This behavior is typical 
of that observed In all the experiments. A slight amount of low frequency 
growth may also be noticed even though the spectral range over which 
comparisons may be made is very limited. At 1 kHz, for example, the 
measured level is about 5 dB above the prediction based on linear theory. 

A comparison of low and high level noise spectra is shown in Fig. 5-2. 

The low level spectra are on the left, the high level on the right. 

The same input noise, a 1/1 octave band centered at 4 kHz (first row in 
Fig* 5-2a) , was used for both sets of measurements, but the input level 
was raised approximately 20 dB for the high level set. 
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FIGURE 5-2a 

COMPARISON OF MEASURED LOW AND HIGH LEVEL NOISE 
SPECTRA (EXPERIMENT AEM4). THE MAXIMUM AMBIENT 
NOISE LEVEL FOR FREQUENCIES GREATER THAN 1 kHz 
IS APPROXIMATELY 25 dB (100 Hz BAND) re 20 /iPa 
R 0 ( 10 kHz) = 7.6 m 

LOW LEVEL: OASPL = 123 dB at 1 m 

HIGH LEVEL: OASPL = 142 dB at 1 m 
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FIGURE 5-2b 

COMPARISON OF MEASURED LOW AND HIGH LEVEL NOISE 
SPECTRA (EXPERIMENT AEM4). THE MAXIMUM AMBIENT 
NOISE LEVEL FOR FREQUENCIES GREATER THAN I kHz 
IS APPROXIMATELY 25 dB (100 Hz BAND) re 20 fiPa 
R o (10 kHz) = 7.6 m 

LOW LEVEL: OASPL = 123 dB at 1 m 
HIGH LEVEL. OASPL = 142 dB at 1 m 












The second row shows the acoustic spectra inside one of the horns (r = -*0.23 m) ; 

the OASPL there was 126 dB for the low level noise and about 145 dB for 

the high level noise, inspection of the remaining spectra in Fig. 5-2 

shows the following. Even for the low level noise, nonlinear propagation 

distortion was not negligible. There is a clearcut growth of the second 

harmonic band and even a development of a third harmonic band. Inside 

the horns the second harmonic band is about 45 dB below the fundamental 

band; at 70 m the two bands are only 20-25 dB apart. The high level 

data shows a much greater effect of nonlinearity. There is a robust 

growth of the entire spectrum above the original 4 kHz band. Even at 

20 kHz, where atmospheric absorption is strong, the growth is impressive. 

Inside the horns the signal at 20 kHz is in the instrumentation noise, 
at least 60 dB below the fundamental band. At 70 m, however, the 20 kHz 
signal is within 30 dB of the fundamental band. To put it another way, 
the 20 kHz level is about the same at 70 m as it seems to be at 0 m. 

Again a small but noticeable low frequency growth may be discerned. At 
70 m the 1 kHz level is closer to the spectral peak in the measurement 
made at high level than in the one at low level. 

Thus far we have said that both diffraction and nonlinear effects 
are responsible for forming the broadband spectrum observed at the end of 
the array nearfield. To quantitatively assess the importance of 
each of these effects, we have computed farfield spectra based on linear 
diffraction theory (see Eq. 3-3). The results of the computation are 
compared with data in Fig. 5-3. (The data are simply a repeat of the last 
three rows In Fig. 5~2b plus a measurement at 18.2 m that was omitted from 
Fig. 5-2.) In both the low and high level cases, the prediction matches the 
data well within the original octave band. At higher frequencies, however, 
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SPL (100 Hz BAND) - dB re 20 /z Pa 











the prediction falls well below the data. On the basis of this observa- 
tion, we conclude that the observed spectral growth in the near field far 
exceeds that which would be due to diffraction alone. Furthermore, because 
the spectral growth is much more pronounced for the high level data, we 
conclude that it is caused by nonlinear propagation distortion. 

Another view of the same experiment (high level measurements only) 
is presented in Fig. 5-4. Here propagation curves for specific 100 Hz 
bands are shown. The data were obtained from the high level spectrograms 
by reading the band levels as a function of range. The solid curves are 
predictions based on linear theory (spherical spreading and atmospheric 
attenuation). The 3.8 kHz band (approximately the spectral peak) decays 
in accordance with the prediction. For the higher frequency bands, however, 
the agreement becomes progressively poorer with both frequency and distance. 
At 70 m, for example, the bands at 8, 14, and 20 kHz are approximately 3, 9, 
and 15 dB, respectively, above the linear theory predictions. It is also 
significant that the divergence between data and linear theory shows no 
sign of slackening, even at the greatest range. In other words, there 
is no evidence in this experiment that small-signal theory will eventually 
triumph. 

Figure 5-5 shows noise spectra recorded to 50 kHz. The input signal 
to the array was a 1/1 octave band of noise centered at 4 kHz. Once more, 
linear theory provides a good description of the behavior of the original 
octave band but a poor one of the behavior of the high frequency noise. 
Moreover, the discrepancy between linear theory and measurement again gets 
progressively worse with both increasing frequency and distance. For 
example, holding distance fixed at 79.3 m, we find the discrepancy is about 
4 dB at 10 kHz, 15 dB at 20 kHz, 30 dB at 30 kHz, and 40 dB at 35 kHz. 
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FIGURE 5-5 

NOISE SPECTRUM AS A FUNCTION OF RANGE 
SEVEN DRIVERS (EXPERIMENT JBL5) 

R 0 (20 kHz) = 7.42 m, OASPL a 140 dB re 20/iPa at 1 







At fixed frequency, say 25 kHz, the discrepancy is about 5 dB at 18.3 m, 

10 dB at 36.6 m, and 20 dB at 79.3 m. In other words, even at the greatest 
distance, the divergence between linear theory and measured data continues 
to increase. 

Figure 5-6 shows results for the first experiment done. A single 
horn driver was used. Because of the broad initial spectrum, the low 
source level, and the low upper limit of the spectrogram, only minor non- 
linear distortion is visible, and it is restricted to the very upper limit 
of the spectrum, 7.5 to 10 kHz. This result is not surprising. Because the 
source level is low, the distortion components are weak. Moreover, the 
broad initial spectrum serves to mask the distortion components produced. 

The good agreement of the linear theory prediction and data within the 
original spectral band (approximately 1 to 7 kHz) shows that the linear 
theory model is suitable for low level signals. 

In Figs. 5-7, 5-8, 5-11, and 5-9 the initial spectral bands are 
approximately the same but the value of source OASPL gradually Increases : 

127 dB, 132 dB, 141 dB, and 144 dB, respectively. As the source level 
is Increased, it is apparent that the distortion occurring in the trans- 
mitter near field is increased. For example, at 9.14 m the 20 kHz level is 
about 35 dB below the spectral peak in Fig. 5-7, but only 19 dB below 
in Fig. 5-9. An observation that is puzzling at first glance is that 
linear theory seems to work better for the higher level experiments than 
for the lower level ones. Compare, for example, the difference between the 
linear theory predictions and the measurements at 20 kHz for the 77 m and 
79.3 m data of Figs. 5-7 and 5-9, respectively. The difference is about 
17 dB for Fig. 5-7 but only 10 dB for Fig. 5-9. The reason for this 
apparent anomaly is that the nearfield distortion, which is much greater 
in the high level experiment, tends to partially mask the distortion 



components produced in the far field. Although the farfield distortion 
components are much weaker in the low level experiment (Fig. 5-7), they 
have comparatively good visibility because the spectrum emerging from the 
nearfield (1.07 m) is much "cleaner'’ than in the high level experiment. 

The results of a narrowband experiment (a 1/3 octave band at 6.3 kHz) 
are shown in Fig. 5-10. The successive major spectral peaks are those of the 
fundamental, second, and third harmonics. Notice how well developed are 
the higher harmonic bands when the noise emerges from the near field. The 
discrepancy between linear theory and measurement is not large in this 
experiment, but it does increase with harmonic number; at 73.2 m it is 
about 4 dB for the second harmonic band and about 10 dB for the third 
harmonic band. The behavior of this discrepancy is in line with predictions 
based on the perturbation solution for tones. A weak difference frequency 
hand centered at about 700 Hz may also be seen. At r = 9.14 m this signal 
was at least partly due to direct radiation from the electroacoustic 
source. The signal is, however, augmented (relative to the linear theory 
prediction) as it propagates from 9.14 m to 73.2 m. 

The data shown in Fig. 5-12a are similar to those in Fig. 5-2 except 
that in Fig. 5-12a the difference in 0ASPL between the low and high level 
experiments is 10 dB, not 20 dB. In addition, the initial bandwidth 
of the noise was somewhat narrower than in the other experiment. The 
spectra in the bottom row of Fig. 5-12a are overlays of spectra measured 
inside the horn and at 9.04 m; they show how extensive the nearfield 
distortion becomes when the level is raised just 10 dB. Figure 5-12b 
shows the farfield spectra but only for the high level expeirment. 

The spectra in the right column were obtained by applying linear horn 
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FIGURE 5-6 

NOISE SPECTRUM AS A FUNCTION OF RANGE r 
SINGLE DRIVER (EXPERIMENT JBL1) 

Rq (10 kHz) = 0.53 m, OASPL = 121 dB re 20 /i.Pa ot 1 I 
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FIGURE 5-7 

NOISE SPECTRUM AS A FUNCTION OF RANGE r 
SINGLE DRIVER (EXPERIMENT JBL2) 

R 0 (20 kHz) = 1.06 m, OASPL = T27 dB re 20 pPa at 1 m 
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FIGURE 5-8a 

NOISE SPECTRUM AS A FUNCTION OF RANGE r 
THREE DRIVERS (EXPERIMENT JBL3) 

R (20 kHz) = 3.18 m, OASPL = 132 dB re 20 /iPa at 1 m 
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FIGURE 5-9a 

NOISE SPECTRUM AS A FUNCTION OF RANGE r 
SEVEN DRIVERS (EXPERIMENT JBL4) 

R 0 (20 kHz) = 7.42 m, OASPL = 144 dB re 20 fjL Pa at 1 m 
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FIGURE 5-9b 

NOISE SPECTRUM AS A FUNCTION OF RANGE r 
SEVEN DRIVERS (EXPERIMENT JBL4) 

R 0 (20 kHz) = 7.42 m, OASPL » 144 dB re 20 /iPa ot 1 m 
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FIGURE 5-l2a 

COMPARISON OF MEASURED LOW AND HIGH 
LEVEL NOISE SPECTRA (EXPERIMENT AEM3) 

FOR THE DATA IS THE LAST ROW THE r - -0.23 m SPECTRA 
HAVE BEEN SHIFTED UPWARD FOR COMPARISON PURPOSES 

R 0 (10 kHz) s 7.6 m 

LOW LEVEL: OASPL = 132 dB at 1 m 

HIGH LEVEL: OASPL - M2 dB at 1 m 
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and diffraction theory to the spectrum measured inside the horn 
(r » -0.23 m). The 14-20 kHz portions of the computed spectra are 
too high because they are based on nonacoustic signals in the 

starting spectrum. That is, the 14-20 kHz portion of the spectrogram at 
-0.23 m represents instrumentation noise, not sound. Even discounting the 
erroneous prediction at the high end of the spectrum, however, comparison 
of the measured and predicted spectra at 9.04 ra shows that much distortion 
occurred in the near field. In the farfield the discrepancy continues to 
grow. For example, the discrepancy at 10 kHz is successively 14 dB, 17 dB, 

19 dB, and 21 dB as range increases. No sign of any leveling off to a 
constant discrepancy appears. 

The experiment depicted in Fig. 5-13 is similar to that in Fig. 

5-2, but at an intermediate source level. The stairstep appearance of the 
nearfield spectra is due to a relatively good resolution of the harmonic bands. 
The resolution breaks down as distance increases because the buildup of inter- 
modulation components fills in most of the sharp breaks in the spectrum. 
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CHAPTER 6 


DISCUSSION AND ANALYSIS OF EXPERIMENTAL RESULTS 

The experimental results are discussed and interpreted in this chapter, 
and some attempt is made to provide theoretical explanations for the ob- 
servations. The discussion in Section A is general and covers both near- 
field and farfield. In Section B the nonlinear theory (for the farfield) 
described in Chapter 3 is tested by comparing predictions with measurements 
from Chapter 5. 

A. DISCUSSION OF EXPERIMENTAL RESULTS 

1. Nearfield Measurements 

Four observations may be made regarding the nearfield measurements. 
First, a tremendous growth of the high frequency end of the noise spectrum 
occurred in the array nearfield. See, for example, Fig. 5-12. Second, 
calculations showed that the growth of the high frequency end of the 
spectrum far exceeded that expected on the basis of linear diffraction 
theory. See Figs. 5-3 and 5-12. Third, on the other hand, the most 
intense portion of the spectrum was well described by the linear diffraction 
theory calculation. This observation is not inconsistent with the previous 
two. Only a small decrease, e.g., a small fraction of a decibel, in the 
intense portion of the spectrum is required to produce significant high 
frequency growth. Fourth, the high level data showed much more high 
frequency growth than did the low level data (see, for example. Fig. 5-2). 

On the basis of these observations, we conclude that the high frequency 
growth is a result of nonlinear propagation distortion. The fact that the 
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spectrum becomes so broad so rapidly is not surprising. Although the pro- 
pagation distance is short, typically 9 m or less, the amplitude remains 
high while the wave is in the near field. 

2. Farfield Measurements 

The nonlinear effects noticed in the farfield measurements were 
similar to those observed for the nearfield data. In all cases high 
frequency growth was evident; in some cases It was strong. Linear 
theory worked well for the intense, middle portion of the spectrum. 

A very limited amount of low frequency growth took place. The distortion 
that occurred was therefore of the Stage I type. The noise in our 
experiments was not intense enough to produce Stage II distortion. 

As noted in Chapter 1, it is commonly believed that while nonlinear 
effects may be important close to a jet engine, where the noise is very 
intense, at some large distance they become unimportant because the SPL 
has reached the "small-signal range." Linear theory may be relied on 
thereafter. This premise is now tested, insofar as it can be tested 
by our experimental data. Consider, for example, Fig. 5-11, which shows 
a discrepancy of about 15 dB between the 20 kHz levels measured and pre- 
dicted at r *= 79.3 m. The linear theory prediction is of course based on 
using 9.14 m as the value for r^, i.e., taking the spectrum at 9.14 m as the 
starting spectrum. If r Q is doubled, that is, if the linear theory pre- 
diction is begun with the spectrum at 18.3 m, one finds that the discrepancy 

at r ■ 79.3 m is reduced only a little, to 12 dB. Even if r is doubled 

o 

again, which means starting the prediction with the 36.6 m spectrum, the 
discrepancy at r « 79.3 m is still 9 dB. In other words, at no distance 
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within our experimental range did divergence from linear theory predictions 
cease. In Chapter 5, the data shown in Fig. 5-4 were also used to make this 
point . 


Although our experimental evidence supports the notion about the high 

•j^ 

frequency noise that, in Pernet's words, "once nonlinear, always nonlinear," 
we have no measurement beyond approximately 80 m. It might be argued, 
therefore, that if we had made measurements at 200 m, or 500 m, or 1000 m, 
we would eventually have encountered small-signal behavior. There is, 
however, sound theoretical basis for expecting that small-signal behavior 
will never be established, no matter how great the distance from the source. 
The theoretical reasoning is similar to that in Chapter 2 except that here 
no specific frequency dependence for a is assumed. Consider sinusoidal 
spherical waves whose amplitude at the reference distance r^ is P-^q* If 

the source level is not too high, the second harmonic pressure amplitude 

5 27 

p^ is accurately given by the perturbation solution * 
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where a^ and a ^ are the small-signal attenuation coefficients at the 

fundamental and second harmonic frequencies, respectively. The asymptotic 

form of Eq. 6-1, valid for (&2 ~ 2a_^)r >> 1, is 
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Equation 6-2 is remarkable in two respects. First, the spreading factor 
-2 -1 

seems to be r , not r . Second, the apparent absorption coefficient is 
2a^, not a 2 * Since generally > 2a^ (a 2 = 4a^ for a thermoviscous medium), 
the nonlinearly generated second harmonic eventually decays more slowly than 


* 

D. 


F. Fernet, personal communication to D. T. Blackstock (1977). 
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would have been expected from linear theory. Moreover, no distance is ever 

reached at which the small-signal behavior r^expCa^r) ] is established. 

The third harmonic is found to behave similarly. At great distance its 

-3 -1 

amplitude varies as r exp(-3ct^r), not as r exp(-a^r). (As noted in Chapter 2, 

2 3 

these results are equivalent to ~ p^ and P 3 * w ^ ere P^ 

local amplitude of the fundamental.) Because we have here considered 
sinusoidal signals, not noise, the analysis is only suggestive. At the same 
time there is no reason to expect the high frequency noise components to have 
markedly different properties. 

B. COMPARISON OF DATA WITH NONLINEAR THEORY 

In this section we compare measured spectra with those computed using 
the nonlinear theory described in Chapter 3. In particular, comparison 
Is made with the data from Experiment AEM3 (see Fig. 5-12). The time 
waveform shown in Fig. 6-1 was constructed from the 9.04 m spectrum by 
the method described in Chapter 3. The sampling interval Af was 50 Hz 
so that the period of the constructed waveform, 1/Af, was 20 msec. The 
constructed waveform was then used as the source waveform in the computa- 
tional procedure described in Chapter 3 0 A comparison of measured spectra 
with the computed ones is given in Fig. 6—2; measured spectra are on the 
left, computed ones on the right. The agreement between measured and 
computed spectra is good for the spectral peak. For the higher frequency 
regions, however, the agreement becomes progressively poorer with both 
increasing frequency and range. The poor agreement at the higher frequencies 
is typical of that found when other comparisons with data were made. 

In general, predictions based on the nonlinear model were at best only slightly 
better than those based on linear theory. 
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FIGURE 6-1 

TIME WAVEFORM CONSTRUCTED FROM THE 9.04 m SPECTRUM 
OF EXPERIMENT AEM3 
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FIGURE 6-2 

COMPARISON OF MEASURED AND COMPUTED NOISE SPECTRA (EXPERIMENT AEM3) 

A NONLINEAR THEORY WAS USED TO OBTAIN THE COMPUTED SPECTRA (SEE TEXT) 









We feel that the failure of the predictions based on our model of 
nonlinear theory is not due to any error in nonlinear theory per se . 

After all, direct application of the theory gave predictions that were in 
excellent agreement with measurements in Phase I. It seems likely that 
our model failed because our method of constructing the source waveform 
was too simple. It will be recalled that the source spectrum, e.g., 
the 9.04 m spectrum shown in Fig. 6-2, was represented as a group of 
closely spaced line components, each having a random phase. No other 
property of the source noise was specified. Yet the source noise does 
have a special character by virtue of the distortion it has suffered in 
traveling through the horn and through the nearfield. In particular, the 
spectral region above the original noise band (the original band was 
2-4 kHz for the noise shown in Fig. 6-2) is made up of distortion components 
(harmonics and intermodulation products) generated by nonlinear interactions 
in the original band. There is thus a closer connection between the high 
frequency components and the original band components than there would 
have been had the high frequency components been generated Independently. 

One property that ought to be considered, for example, is the phase 
coherence between components in the original band and those at higher 
frequencies. By neglecting phase coherence when we constructed the source 
waveform, we eliminated the possibility that subsequently generated 
higher frequency components could add in phase to components already 
present. This would help explain why the predicted spectra fall below 
the measured spectra at high frequency (Fig. 6-2). 

The simplest way to correct the error in our prediction model would 

be to use a direct recording of the noise signal at r = r as the input 

o r 

waveform for the computer program. 
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Previous conclusions in this chapter regarding the role of nonlinear 
effects in our experiments are, of course, in no way affected by the success 
or failure of our nonlinear propagation model. The conclusions were estab- 
lished independent of such a model and hence stand alone. 



CHAPTER 7 


APPLICATION TO JET NOISE 


The relevance of our results for actual jet aircraft noise is discussed 
in this chapter. Because the noise in our experiments is approximately 
one decade higher than typical jet noise, some scaling must be done before 
comparisons can be made. 


A. SCALING LAWS FOR NONLINEAR PROPAGATION DISTORTION 

As was pointed out in Chapter 2, when ordinary absorption is not 
important, the distortion of two different sinusoidal waves is the same 
if the amplitude- frequency product is the same for each (see the dis- 
cussion following Eq. 2-1). We now seek a generalization of this rule 
for waves of arbitrary waveform. 

The derivation given here is similar to that in Ref. 13. We start 
with the Earnshaw solution, Eqs. 3-7. For convenience Eq. 3-7b is 
repeated here 

t* = t - Bu(t)x/c 2 . O" 1 ) 

o 

It will be recalled that t' is the retarded arrival time of a wavelet 
whose particle velocity is u. The parameter t is the value of t’ at x=0. 

For simplicity of notation in what follows we use the unprimed symbol 
t in place of t f . In Fig. 7-1 the distortion of a wave is traced. Two 
particular wavelets are identified by their particle velocities ^ and 
u 2 and their values of t (t lQ and t 2Q , respectively). Application of 
Eq. 7-1 to these two wavelets gives 


and 


H - ho " Bu i x/c o 2 


h = ho - Su 2 x/ V 


(7-2) 
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( 7 - 3 ) 


The time separation of the two wavelets At * t^ - t^ is therefore 

At * At - gAux/c 2 j 
o o 

where At Q = t£Q - t^Q is the original separation and Au ■ - u^. The 

fractional change in separation of the two wavelets 6 = (At Q - At)/At Q is 

6 * SAux/c 2 At . 

o o 

The value of <$ tells us how close one wavelet has come to the other and is 
therefore a measure of the distortion of the wave. For instance, 6=1 
implies the second wavelet has caught up to the first one. It will be 
seen that 6 can be doubled by doubling the amplitude of the wave (because 
u is then doubled) or by time compressing the wave by a factor of two 
(because At^ is then halved). Time compressing by a factor of two means 
doubling the frequency of every component in the spectrum. Alternatively, 
if the amplitude is doubled and the frequency halved, the degree of dis- 
tortion will remain the same. Thus for random signals as well as for 
periodic ones, it is the product of amplitude and frequency that determines 
the degree of distortion. For spherical waves the expression for 6 may be 
found by replacing x in Eq. 7-4 by r^lnCr/r^). The frequency-amplitude 
scaling is not, however, affected by this change. 

B. COMPARISON OF MEASUREMENTS WITH JET DATA 

Spectra from two of our propagation experiments, JBL9 and AEM3, are 
compared with spectra of a KC-135A aircraft (ground runup, single engine, 
measurement angle 30* off the jet axis) in Fig. 7-2.* The distance for the 
comparison was 250 ft, or 76.2 m (the data had to be extrapolated backward 
about 3 m, a correction of 0.3 dB, in order to make the comparison). The spectral 

*The KC-135A data was kindly furnished by John N. Cole, Aerospace Medical 
Research Laboratory, W-PAFB, Ohio. 
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SPECTRUM LEVEL - dB re 20 /iPa/H 



FREQUENCY - kHi 


FIGURE 7-2 

COMPARISON OF ACTUAL JET NOISE SPECTRA WITH SCALED SPECTRA 
(SCALED DOWN BY A FACTOR OF 10 IN FREQUENCY, SCALED UP 20 dB 
IN LEVEL) FROM EXPERIMENTS AEM3 AND JBL9 


RANGE = 76.2 m 




peaks for our data were approximately a decade above those for the KC-135A. 
For ease of comparison, therefore, we have scaled our spectra down in 

frequency by a factor of 10 and up in spectrum level by 20 dB, in accordance 
with the scaling law given in Section A. If it is assumed that the two 
noises had similar statistical properties at their respective sources, 
for example, if both were Gaussian, then it is reasonable to compare them 
for their susceptibility to nonlinear propagation distortion. 

It will be seen that the level of KC-135A noise is roughly 10 dB higher 
in the mid- and high-frequency regions than our noise. By demonstration, 
our noise was definitely affected by nonlinearity. The implication is 
that even stronger nonlinear effects were at work during the propagation 
of the KC-135A noise. Moreover, although the KC-135A is a very noisy 
aircraft, many other current aircraft produce noise whose spectrum levels 
are higher than our scaled spectra. One therefore concludes that nonlinear 
effects are probably common in jet noise. 
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CHAPTER 8 


SUMMARY 

A series of experiments has been carried out on the propagation of 
finite-amplitude noise outdoors. The purpose of the study was to determine 
the extent to which nonlinear effects influence the propagation of noise 
from a controlled source in the outdoor environment. A secondary purpose 
was to develop, if possible, a theoretical model for the noise propagation. 

A ground -mounted electroacoustic source transmitted broadband, octave 
band, or 1/3 octave band noise in the frequency range 2-10 kHz, The 
source level (overall sound pressure level) of the noise was in the range 
121-145 dB re 20 yPa at 1 m. The propagation path was vertical and parallel 
to an 85 m tower, whose elevator carried the traveling microphone. The 
maximum propagation distance was about 80 m. The experiments were done at 
night during the months June through September 1977. The meteorological 
conditions (ground level) were as follows: temperature range 23-31 °C, 

relative humidity range 55-90%, and wind speed range 0-24 km/h. 

The measurements were compared with theoretical predictions based on 
linear theory. Spherical spreading, atmospheric attenuation, and, as 
appropriate, source diffraction were accounted for in these predictions. 
Predictions based on nonlinear theory were also attempted, but the partic- 
ular model was based on an assumption about the noise distortion in the 
transmitter nearfield. As it turned out, this assumption was not justified 
in our experiments. 

Use of an amplitude-frequency scaling law made it possible to compare 
the noise from our experiments with noise from a KC-135A aircraft. The 
importance of nonlinear effects in actual jet noise could thus be estimated. 
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The primary conclusions are these. 

1. A strong generation of high frequency noise caused by nonlinear 
effects was found in all the high intensity noise experiments. A 
very limited amount of low frequency noise was also generated. The 
intense, middle part of the spectrum deviated little if any from 
expectations based on linear theory. These observations indicate 
that although shocks formed in the noise waveform, distortion did not 
reach the stage at which shock merging was important. 

2. The spectral distortion occurred in both the transmitter nearfield 
and far field. Moreover, the distortion in the near field was over and 
above the spectral changes associated with diffraction. 

3. At no measurement point was small-signal behavior established 
for the high frequency noise. Theoretical calculations for tone 
signals support the proposition that the nonlinearly generated high 
frequency noise never achieves a farfield where small-signal behavior 
is established. 

4. Comparison of the measured spectra with predictions from a 
model based on nonlinear theory showed poor agreement for the high 
frequency noise. The failure of our predictions was not due to an error 

in the nonlinear theory per se, but rather due to an inadequate description 
of the source noise waveform. In the future, use of a directly recorded 
input waveform should be tried. 

5. Comparison of our scaled experimental measurements with actual 
jet spectra show that the spectrum levels encountered in our ex- 
periments are well within the jet noise range. Indeed, the noise 
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measured from a KC-135A jet is roughly 10 dB higher in spectrum 
level than our scaled noise. One therefore concludes that nonlinear 
effects are probably common in jet noise. 

The report also contains some information about outdoor propagation 
of finite-amplitude tones. Tone experiments done as a forerunner to those 
on noise are reviewed. Additional data on tones taken during the course 
of the noise study are also described. Theoretical analyses of propa- 
gation of very weak and very strong tones are presented. 



APPENDIX A 


THEORETICAL RESULTS FOR FINITE-AMPLITUDE 
SPHERICAL WAVES 2 

I. WEAK WAVES - A PERTURBATION SOLUTION OF BURGERS' EQUATION 


Our purpose here is to obtain a theoretical prediction valid for 

weak spherical waves in a thermoviscous fluid. Several other researchers 

5 27-29 

have obtained results valid for the weak-wave problem. ' 

Pernet and Payne, ^ for example, extended the phenomenological approach 

30 29 

Thuras et al had used for plane waves. Blackstock and Willette 

started with the Burgers' equation for spherical waves and obtained a 

perturbation solution valid to third order in the source Mach number e . 

The work we present here is an extension of the perturbation solution to 

fifth order in e. 

The Burgers' equation for spherical waves is (see, for example. 

Ref. 31) 

... .? 

0 „ 3l 


l dOp) a 3 p 
r ~ Sr 


a ° P P r> £E . 

a ) 2 at ' 2 p - c - d ^ 


(A-l) 


o o 


where a is the thermoviscous absorption coefficient at angular frequency U) 

and t ,,a t-(r-r )/c Is the retarded time. Let the boundary condition be 
o o 


p(r ,t) = P in sin cot 


(A- 2) 


10 

If the dimensionless quantities V=(r/r ) (p/p c 2 ), y=wt', and £=3kr 

o o o o 

are introduced, Eqs. A-l and A-2 become 

, 2 . 


dv a 
3? ' T 2 


3‘v t .. dv 
3y 2 " r 


<5 y 


and 


V(r ,y) = € sin y 


(A-3) 


(A-4) 


respectively. We seek a solution of Eq. A-3 in the form of a perturbation 
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series. That is, we suppose a solution of the form 


V - E «V n) 

n=l 


CA-5) 


exists. Substituting Eq. A-5 into Eq. A-3, we obtain 

d 2 v (n) £ y- v (i) ^ 


f n [V n) 

£ 6 \SF- - 
n=l L 


a 


v' 


: ] 


= o 


CA-6) 


dy~ " i+j=n 

In order that Eq. A-6 be satisfied for artitrary values of e, the bracketed 
term must vanish for each value of n. The result is the following infinite 
set of linear. Inhomogeneous differential equations: 


3v 

"3? 


- a » — 


-i E 


v 


(i) 


U) 




(A- 7) 


dy^ x i+3=n 

The equations must be solved in sequence. For example, the source 
(inhomogeneous) terms for the second order solution depend on 

The source terms for the third order solution depend on and 

We now proceed to solve the first five equations. 

The first order equation is 


3v 

“3? 


- a * — = o 


(A-8) 


V 


The solution satisfying the boundary condition Eq. A-4 is 

,,(l) . 

V = e sm y (A-8a) 

This is the familiar solution of linear acoustics. The second order 
equation is 


"3r 


- a 


sV 2) _ a „(i) 3v 


(i) 


V 


‘ 2r 


-2a(r-r ) 


sin 2y 


(A-9) 
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whose solution Is 


r(2) _ L 


-lja(r-r ) 


I 22 sin 2 y > 


(A-9a) 


where 


2a(r 1 -r ) 


■22 


* o 


dr ' 


We shall have occasion to write down several integrals of the form 

y*F(r)r ^exp(ar)dr and adopt the following convention. The integral 

th 

I is associated with the contribution of the m order solution to the 
mn 

n^ harmonic. 

The third order equation is 

acv (1 V 2) ) 


dv 

" 5 ? 


(3) d 2 y (3) . 

- a r — = - 2 - 






4 


-5cc(r-r ) 


I 22 e 


° [3 sin 3y - sin y] . (A-10) 


Its solution is 

>2 -a(r-r ) 


(3) - - f- 


V 

where 

and 


Si 6 


2 -9a(r-r ) 

sin y + I e ° sin 3y > (A-lOa) 


'31 


■/, 

■/, 


r I 22 (r') 6 


-^(r'-r ) 


r I 22^ r '^ 6 


^(r'-r ) 


dr' 


dr' 


We see that the third order solution contains not only a third harmonic 
term but also one that represents a correction to the amplitude of the 
fundamental. This correction denotes loss caused by nonlinear effects. 
The fourth order equations and solution are listed below. 
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dV 4 ) 


- a 


dy c 


= i f v (2) 9 v (2) + 3C 

r ^ <3y dy J 


= g 


-2o:( r-r ) 


•r I V _I 3X e 


- 3 I 33 8 


-10a( r-r 


'•’) ■* 


sin 2y 


*( 


2 -8a(r-r Q ) -10a(r-rJ 

122 ® © 


o) ) si 


sin 4y 


- 


-4a(r-r ) 
(i^o+ji'i.o) 8 0 


h2 rj ^ 42 } e sin 2y 

-l6a(r-r ) 


+ ( x 44 + 611 44) e U sin V 

2a(r 1 -r ) 




-42 


I* 


42 


-44 


-/, 

■/, 

■L 


r I (r' ) e 


dr r 


-6a( r 1 -r ) 

rl 33 (r')e 


2, 8“(r’-r ) 


dr 1 


r I 22 (r ’ ) e 


dr’ 


L 

j r 


r I«(r 1 ) e 

T • _ I 

1 44 


6a(r f -r ) 


dr 1 


* (A-ll) 


(A- 11 a) 


Because of the rather formidable algebra involved, the only part of 
the fifth order solution that was calculated was the contribution to the 
fundamental. This contribution is denoted . The results are 
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listed below. 



(A- 12) 


(A- 12 a) 


‘51 


I, 


o L 


I k2 (r') + 31 Vo(r' ) + 


42' 


I 22 (r') I ;i (r') 


-^a(r'-r ) 
o' 


dr ' 


' 51 = L 


r I 22 (r’) I 55 (r') e 


-12a(r '-r ) 
o 


dr' 


As may be seen from the higher order results, the solutions are easy to 

write down in terms of the integrals I . It is the calculation of these 

mn 

integrals which is indeed a arduous task. 

Let us now put our results in terms of the normalized harmonic 

amplitudes B =(r/r )(p /p Defining 0 =e£, we obtain 
n on jlu o 


B, 


-a( r-r ) 
N o 
e 


— r-r 

' o 
e 



B 


3 


-9«(r-r o ) 
e L ' 



) 


CA-13) 


CA-14) 


CA-15) 
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and 


CA-16) 


= e 


-l6ct(r-r Q ) 


E 


a 5 
o 

-r 


di^+ei'ht) 


1*4 


'] 


We note here that in the limit as ar+O, the expressions for the 

given in Eqs. A-13 to A-16 reduce to the corresponding expressions 

found using lossless theory directly. The details are given in Ref. 2. 

In many cases where a perturbation solution is useful, some 

distortion is evident in the source waveform V(r ,y). Commonly 

o 

this distortion is primarily second harmonic (especially if it is due to 
nonlinear propagation distortion in the nearfield). We shall now 
determine how the presence of an "initial" second harmonic signal changes 
the lower order perturbation results. Suppose that, in place of Eq. A-4, 
the boundary condition is 


V(r Q ,y) = 


sin y + a sin 


(2y+q>)J 


(A-17) 


where a and <j> are the relative amplitude and phase of the second harmonic, 
respectively. The first order perturbation solution satisfying this 
boundary condition is 


V 


( 1 ) 


-a(r-r Q ) 


sin y + a e 


-4a(r-r ) 


sin(2y+cp) 


(A-18) 


The second order equation, Eq. A-9, has a source term proportional to 

V (1) (3V (1) /3y). Hence, it is easily seen that contains contributions 

( 2 ) 

to the first, second, third, and fourth harmonics. The expression for V 
is 


V 



-a(r-r Q ) 

e 


X 21 Sin(y+Cp ) 


+ 


L 

2 


-4a(r-r Q ) 

e 


I 22 Sln 2y 


+ F(r) sin(3y+<p) + G(r) sin(4y+2cp) , 


(A-19) 
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where 


■i: 


-lja(r *-r ) 
r ' o 
e 


21 / r 

o 


dr’ . 


r(2) 


We note here that the second harmonic component of V is not affected by 
the "initial" second harmonic. To second order, therefore, the expressions 
for p^ and p^ are 


P 1 = p io e 


-a:(r-r ) 
o' 


1 + 




1/2 


X 21 ‘ a °o X 21 cos ^ 


, (A-20) 


and 


P 2 = P 10 e 


-ta(r-r„) / , ff 
z o 


1/2 


a + 


nr X 22 + aa o J 22 COS * 


(A-21) 


A comparison of Eqs. A-20 and A-21 with experimental data is given in 
Chapter 2. Further comparison of Eq. A-21 with measured data is given 
in Appendix D. 

The asymptotic values (ar-+°° ) of the harmonic amplitudes p^ are 
particularly interesting. From Eqs. A-13 to A-16, the asymptotic values 
for the first four harmonics are 


and 


^ p l^asm P 10 \ 


-a (r - r ) 
e 0 

* 

(A- 2 2 a) 

^• P 2^asm P 10 ^ 


r\ -2a (r - 
0 e 

r o ) • 

(A-22b) 

^ p 3^asm P 10 ^ 

4 

L) 2 e- 3oCr - 

r o^ > 

(A-22c) 

^ p 4^asm P 10 ( 

4 

v 3 

r\ -4a(r - 

4 e 

*0) . 

(A-22d) 


As in the plane wave case (see Chapter 1) the apparent absorption coefficient 
2 

is na, not n a. In addition the apparent geometrical spreading is as r , 


not r 


-1 


The harmonic decay rates given in Eqs. A-22 are, however, slower 


than those predicted using linear theory. That is, for n>l, r exp( nar) 
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-1 2 

ultimately represents a slower rate of decay than r exp(-n ar) . In other 
words, the higher harmonic sound never establishes a farfield where 
traditional small-signal laws take over. 

In conclusion, we have presented a fifth order perturbation solution 
of Burgers' equation. Although the integrals involved become increasingly 
complicated for the higher order results, the low order results should 
be useful for a great many weak-wave problems. 


II. STRONG WAVES - SOLUTION OF THE AMPLITUDE DECAY RATE EQUATION 
33 

Rudnick originally proposed a model for the decay of the peak 

amplitude of a plane, sawtooth wave. The total decay rate for the peak 

amplitude was taken to be the decay rate due to ordinary absorption plus 

the decay rate due to nonlinear effects. We call this an "amplitude 

decay model" to distinguish it from various intensity decay models (see, 

for example. Ref. 21). Blackstock developed analogous equations for the 

* 


fundamental pressure amplitude for both plane and spherical waves, 


It was assumed that the decay rate dp^/dr of the fundamental pressure 
amplitude p^ is the sum of the decay rate due to absorption (dp^/dr^^g “ -a P^ 
and the decay rate (dp., /dr)- of a spherically spreading sound wave of 
finite amplitude. The latter decay rate may be found as follows. If the 
boundary condition is given by Eq. A-2, the value of p^ in the sawtooth 

32 

region (ct >3 ) in the absence of absorption is (from the weak-shock solution ) 

r 2p. 


= TO 

^1 r 1 + cr 


(A-2 3) 


From this equation we obtain 


D. T. Blackstock, unpublished research (1971). 
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(S') 


r Q 2p 10 P€ta- o 2p 10 


1 + a 

f .a. r 


r 2 ( 1+cr ) 2 


1 

= - - P 


£ek 


r "I 2p l0 


By supposition the total decay rate dp^/dr is then 


te P] 

dx 


d P] 

dr“ 


abs 


'dp. 

dr 


f . a. 


(A-24) 


^ 1 pek 2 

- 7 p i - p - 


l 


(A -25) 


We now wish to solve Eq. A-25 subject to the condition that the solution 
reduce to the weak-shock solution (Eq. A-23) as ct^O. The form of the 
solution for the plane wave case suggests that we try a solution of the form 

-a( r-r ) 


o A e 

P 1 = r 1 + f ( r ) 


(A-26) 


where A is a constant and the function f(r) is to be determined. Sub- 
stitution of Eq. A-26 into Eq. A-25 leads to the following form of f(r): 

A £3ekr Ctr 


f(r) - - 


2p 


o o 
e 


10 


/ ’ca 

or 


ar 

e 


d\ 


(A-27) 


The choice A=2p^^ yields the weak-shock solution, Eq. A-23, as a-^0. 

If the definition of the exponential integral E^(a) (s ee > f° r example. 
Ref. 34), 


E 1 (a) 


f°° 

J a * 


d K 
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is used, the expression for p^ becomes 

-Cd(r-r ) 
_ o 


o 


^10 e 


1 + £€kr o e °[E^(o!r o )-E 1 (ar) 1 

Equation A-28 is simple to evaluate in practice. The exponential 


integral E^(ar) may be read from tables or evaluated by certain series 
representations. See, for example, Ref. 34, which gives a six term 


approximat ion 

5 

E 1 (ar) = -ln(ar) + £ a^ar) , (A-29) 

having a maximum percentage error of 9x10 ^ for the range O^ar^l. 

Equation A-28 was tested by comparing it with freefield propa- 
gation data taken using a siren operating at 6.1 kHz in air (see Ref. 1 
or Chapter 2 for a description of the siren). Because of the relatively 
high acoustic power output of the siren (approximately 600 W of acoustic 
power with a 40° beamwidth between the 3 dB down points) , the characteristic 
sawtooth waveshape was already in evidence at a propagation distance of 
approximately 2 m. Our measured boundary condition was therefore a 

sawtooth wave at a given distance r from the source. To use Eq. A-28 

m 

we must compute from our measured data the amplitude p^Q and effective 

source radius r of a sinusoid which would yield the measured sawtooth 
o 

wave at range r . The weak-shock solution may be used for this purpose 
m 

provided ar<<l. The smallest value of r for which Eq. A-23 is valid is 

the one for which a=3. Let p. be the measured value of the fundamental 

lm 

amplitude at range r . Combining Eqs. A— 23 and 2—1 (with o=3, r=r , and 

m m 
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p =p, ), we obtain 
1 xm 


and 


r = r exp (-3/ 2 Be kr ) 
o m mm 


10 


= 2p_ r /r 
lm m o 


(A-30) 


(A-31) 


where e = p_ /p c 
m lm o o 

Figure A-l shows the results of a comparison of Eq. A-28 with the 

measured data. The absorption coefficient a was measured by means of 

a small-signal experiment done under similar conditions and found to be 

0.0059 Np/m, a value consistent with the value computed according to Ref. 22. 

3 

The values of p. and r were 2.29x10 ubar and 2.19 m, respectively. The values of 
lm m 


P^q and computed from Eqs. A-30 and A-31 were, respectively, 1.1x10 ybar and 
0.09 m. The solid curve is the linear theory prediciton; the dashed 
curve represents the solution of the amplitude decay equation. The 
latter is seen to provide a good fit to the data. 

In summary the model proposed here, although ad hoc, does seem to 
provide a good description of the experimental data. Further experimental 
verification of this model may be found in Ref. 1 or Chapter 2. The 
model should be useful in spherical wave problems whenever the source is 
strong enough to produce a sawtooth wave. 
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RECEIVED SPL - dB re 20/i Pa 



FIGURE A -1 

COMPARISON OF THE AMPLITUDE DECAY 
RATE MODEL WITH PROPAGATION DATA 
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APPENDIX B 


AMPLITUDE DENSITY OF A FINITE-AMPLITUDE WAVE* 


Theoretical results are given here for the amplitude density of a 
wave subject to nonlinear propagation distoriton. Rudenko and Chirkin 
found an expression for the amplitude distribution of an initially narrow- 
band Gaussian process. ^ They couch the problem in terms of a complete 
statistical description of the noise. We give results here valid for an 
arbitrary waveform prior to shock formation. 

Figure B-l shows the waveform of an arbitrary signal at x=0. The 
probability P(u, u + Au) of observing wavelets having particle velocity 
amplitudes between u and u + Au is 


N 

P(u, u + Au) = ^ At^/T 
i=l 


(B-l) 


where T is the sample length, N is the number of intervals At^ and the 
{At^} are as shown in Fig. B-l. The probability density p(u) is defined by 


p(u) - J im 
Au-+0 


p(u, u + Au) 
Au 


(B-2) 


Substitution of Eq. B-2 into Eq. B-2 leads to the expression 


p(u, 0) = 


1 

T (0) 


2 

i«l 


(- 1 ) 


i+1 


(3u/3t ) 


(B-3) 


i+1 


where the zero in the argument denotes x=0. The factor (-1) is 
necessary to convert the slopes (3u/3t)^, which- are alternately positive 
and negative, into all positive numbers (all the At^ in Eq. B-l are positive) 


*The text of this appendix has been submitted to the Journal of the 
Acoustical Society of America for publication as a Letter. 
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The + sign is to be used if the first term is positive, the - sign if the 
first term is negative. Equation B-3 shows that the initial amplitude 
density is determined by the collection of slopes (3u/3t)^ of wavelets 
whose amplitude is u. Special note must be taken of what is meant by 
sample length. We take T to be the time interval between the two end 
wavelets in the sample, those at t=0 and t*=T(0) in Fig. B-l. The dis- 
advantage of this definition is that the time interval T may shorten or 
lengthen slightly with distance (as the two end wavelets migrate from their 
initial positions in the waveform) . The advantage is that the particular 
wavelets included in the sample are always the same. If T were chosen 
to be a fixed time interval, some wavelets would in general be lost, or 
new wavelets would enter, at the ends. 

Now let the signal travel to point x, distorting as it goes. The 
new density p(u,x) is 

N 

p(u * x) - tF) I * (iu/at' T Y ’ (B ~ 4) 

i=l 1 

where T(x) is the new duration of the sample waveform, t' = t - x/c 
is the retarded time, and c q is the small-signal sound speed. Because of 

the distortion, the wavelets u are at new positions on the retarded time 

axis. Thus in general both the waveform slopes and the sample duration are 

different. 

In order to determine the amplitude density of the distorted wave 
at point x, a description of the waveform there must be found. Given the 
initial waveform 

u(0,t) * g(t) , (B-5) 
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the solution appropriate for plane progressive waves in a lossless fluid is 

u(x,t) » g(r) , (B-6a) 

where 

t « t'+$xu(T>/c 2 , (B-6b) 

o 

and 0 is the coefficient of nonlinearity. Applying Eqs. B-6 

to the waveform of Fig. B-l, we see that the duration T(x) of the distorted 

sample waveform is 


T (x) = T q - <5T , (B-7a) 

where T q has been used as short for T(0), 6T is given by 

6T = 0x[u(T ) - u(0)]/c 2 , (B-7b) 

o o 

and u(0) and u(T q ) are the particle velocity values at the endpoints of 
the initial waveform. The waveform slopes at point x may also be determined 
from Eqs. B-6. One obtains 

1 _ 1 6x 

3u/3t’ 3u/3t “ IP 2 * 

o 

or, because the slope in the t system is the same as the initial slope 
(see Eqs. B-5 and B-6a) 


1 

3u/3t * 


1 

3u(0,t)/3t 



The expression for the density p(u,x) may now be evaluated. 


Substitution of Eqs. B-7 and B-8 in Eq. B-4 leads to 


p(u,x) 



y . (-d 1+1 + ex[(-D w - 1] 

[3u(0,t)/3t] i " 2 (T q - T)c o 2 


(B-8) 


(B-9) 


2 

because the constant -0x/c in Eq. B-8 cancels out for every pair of terms 

o 

in the series. Only if N is odd is there an odd constant left 
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over, and its presence is accounted for by the final term on the right- 
hand side of Eq. B-9. 

Comparing Eqs. B-3 and B-9, we see that 


p(u, x) 


= p(u, 0) 

1 - 6T/T " 
o 


Sxj(-1) N - 1] 

2 (T - 6T)c 2 
o o 


CB-10) 


It may be seen that for a large class of cases p (u,x)=p(u,0) , that is, 

the amplitude density remains constant as the signal propagates. For 

example, if u(0)=u(T ), then 6T vanishes and there is no leftover term 

o 

(because the number N is even). Even if 6T has a nonzero value, p(u,x) 
may be made arbitrarily close to p(u,0) by choosing T q sufficiently large. 

In either case the amplitude density does not change with distance, so 
long as shocks do not form. For example, the sample noise waveforms 
considered in this report are periodic and continuous (see Chapter 3) . 

The amplitude densities for these signals would not, therefore, be 
expected to change with distance. 

A geometrical explanation of the constancy of the amplitude density 
may be found by examining Fig. B-2. The initial waveform is shown 
in Fig. B-2a. Figure B-2b shows the distorted waveform at point x. 

The size of Au compared to u has been exaggerated for the sake of illustration. 
For a fixed value of u and Au, the probability that the particle velocity 
lies between u and u + Au is proportional to the sum of the intervals 
At^ shown in Fig. B-2. As the wave distorts, the interval A^q is shortened 
by an amount porportional to Au. The interval At 20* ^ owever » *- s lengthened 
by the same amount, so that At^Q + At^Q = At^ + At£» Thus the shortening 
of intervals on the left side of waveform peaks is compensated for by the 
lengthening of intervals or the right side of the peaks. 
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This illustration may be extended to show why p(u,x) must change 
after shocks form. The limiting value of At^ is zero, the value reached 
when wavelets in the window on the left-hand side of the peak form a shock. 
The companion interval At^ is not correspondingly limited; it continues 
to increase after At^=0. The sum At^+At£ therefore departs from constancy 
as soon as the wavelets included in At^ form a shock. 

The conclusion that the amplitude density does not change prior to 
shock formation would seem to be at odds with the results of Rudenko and 
Chirkin, who find that the amplitude density of an originally narrowband 

Gaussian noise changes over arbitrarily small distances. The discrepancy 
is only apparent, however, not real. Rudenko and Chirkin deal with 

signals having arbitrarily large peak values (hence arbitrarily large 
slopes) so that shocks form immediately. Our method should be applicable 
to ’’clipped" signals with finite slopes, where the shock formation distance 
has a nonzero value. 
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APPENDIX C 


MEASURED RANDOM FLUCTUATIONS IN SHORT-TERM SPL FOR TONE SIGNALS 

As was mentioned in Chapter 4, considerable averaging was necessary 
to obtain the data presented in this report. Gusting winds and other 
Inhomogeneous medium effects caused severe fluctuations in short-term 
sound pressure level. [By "short-term sound pressure level" we mean that 
the measurement averaging time is large compared to the period of the 
transmitted signal but small compared to the characteristic time of fluc- 
tuations caused by inhomogeneous medium effects.] We present here some recordings 
of the fluctuation In the short-term SPL as a function of time for several 
ranges. The experimental arrangement is the same as that used for the 
noise experiments. The input to the array for these measurements was, 
however, a sine wave of fixed frequency and amplitude. The amplitude was 
low enough that nonlinear effects were not important. 

A typical set of measurements is shown in Fig. C-l. The ordinate is 
short-term SPL in a 50 Hz band (for all measurements reported here the 
filter was centered at the source frequency) and the abscissa is time. 

The quantity 20 is the half-power (3 dB down) beamwldth of the source at 

Hit 

the given source frequency. At 1.05 m the received level is nearly 
constant. At the two greater distances, however, there are significant 
fluctuations in level. Even though conditions were calm, peak to peak (p-p) 
fluctuations of 4 dB exist at 21. '2 m. Figure C-2 shows data taken under 
windier conditions. Again the fluctuations are seen to increase with 
distance, becoming approximately 10 dB (p-p) at 18.3 m. Finally, measure- 
ments taken with a narrower beamwidth source are shown in Fig. C-3. In 
Fig. C-3a a plot of the long-term (30 sec) SPL versus range is shown. The 
solid curve is a prediction based on spherical spreading and atmospheric 
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1.05 



21.2 


60 90 

TIME - sec 


FIGURE C-1 

MEASURED SHORT-TERM (87 msec) AVERAGE SPL AT VARIOUS DISTANCES 
FROM AN 8.111 kHz SOURCE (2 0 HP = 16°). THE ARROWS INDICATE LEVELS 
PREDICTED FROM SPHERICAL SPREADING AND ATMOSPHERIC ABSORPTION. 

CALM CONDITIONS 
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SPL (50 Hz BAND) - dB re 20 ftPa 



TIME - s 


FIGURE C-2 

MEASURED SHORT-TERM (87 msec) AVERAGE SPL AT VARIOUS 
DISTANCES FROM A 6.45 kHz SOURCE (2 0 HP = 21°). 

THE ARROWS INDICATE 30 sec AVERAGE LEVELS. 

WINDS OUSTING UP TO 20 km/h. 



SPL (50 Hz BAND) - dB rt 20 M Po 



RANGE - m 
(a) 


RANGE - m 



TIME - sec 


(b) 


FIGURE C-3 

(o) LONG-TERM AVERAGE (30 sec) SPL AS A FUNCTION OF RANGE, 
(b) MEASURED SHORT- TERM (87 msec) AVERAGE SPL AT VARIOUS 
DISTANCES FROM A 3.SS kHz SOURCE (2 0 HP = 8.8°). 
WINDS GUSTING TO 25 km/h 
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attenuation. In Fig. C-3b the fluctuations In level are shown at several 
ranges. Despite the fact that the fluctuations in level increase to 
approximately 17 dB (p-p) at the greatest distance, the 30 sec average 
levels (Fig. C-3a) are in excellent agreement with the linear theory 
curve. This observation is the basis for comparing our measured spectra 
with theoretical predictions for a homogeneous medium (see Chapter 5) . 

In other words, if sufficient averaging times are used, a homogeneous 
medium theory seems to provide valid theoretical predictions for our 
experiments. 
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APPENDIX D 


PROPAGATION EXPERIMENTS WITH TONES 

In the course of evaluating the AEM array for use as a noise source, 
we made some propagation measurements with tone signals. The experimental 
arrangement was basically the same as that described in Chapter 4. The 
results given here were previously reported in Ref. 35. 

The AEM array is a rectangular array composed of 20 individual square 
exponential horns, each driven by a separate 35 W horn driver. Each 
horn flares exponentially from a round throat 2.5 cm in diameter to a 
square mouth of side length 11.4 cm. The horn length is 35 cm. The 
horns are bolted together to form a 4 element by 5 element array with sides 
of 45.7 and 57.2 cm, respectively. The array is designed to be used 
as a fog horn. Its frequency response is approximately 250 Hz to 5400 Hz. 

The results of six propagation experiments (AEM 5, 6, 1, and 2 in 
Table 4.1) for which the electrical input was a pure tone are shown 
in Figs. D-l to D-6. The experiments were performed at three frequencies — 
1.5, 3.55, and 5.33 kHz — each corresponding to a relative peak in the 
array response. For each frequency, measurements were made at two different 
source levels,, one 16-19 dB lower than the other. In each case the 
harmonic distortion components, as well as the fundamental, were measured 
as a function of propagation distance. The fundamental is identified by 
the symbol p^, the second harmonic by an< * so on. Each solid curve 
is a linear theory prediction (spherical spreading and ordinary absorp- 
tion). The effective source radius r was taken to be the 

o 

first measurement distance greater than the Rayleigh distance R^f^), 
where f^ is the transmitting frequency. See Table 4.1 for the actual 
values of R Q (f^) for each experiment. For values of the range r less than 
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FIGURE D-l 

PROPAGATION CURVES FOR AEM ARRAY 

AEM5 (LOW LEVEL) 
f, = 1.5 kHz 
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FIGURE D-6 

PROPAGATION CURVES FOR AEM ARRAY 

AEM2 (HIGH LEVEL) 
f, = 5.33 kHz 




R q , the solid curve is an extrapolation of the far field data. The dashed 
curves are predictions of the second harmonic amplitude based on a 
perturbation solution of Burgers' equation (Eq. A-21). The phase angle (f> 
between the fundamental and second harmonic was taken to be zero. [The 
phase angle for the data in Fig. D-6 was measured and found to be approximately 
zero. Time waveforms were not measured for the data of Figs. D-l to D-5.] 

Three conclusions may be drawn from the data in these figures. First, 
for the low level tones, harmonic distortion is small. For the high level 
tones, however, many harmonics are generated along the propagation path 
and some are prominent. For example, in the high level test at 5.33 kHz, 
the second harmonic comes as close as 9 dB to the fundamental. Second, 
in all cases the measured value of p^ was in good agreement with the linear 
theory prediction. These two results are not inconsistent. Only a small 
decrease in SPL of the fundamental (usually only a fraction of a decibel) is 
required to produce significant harmonic structure. Third, the second 
harmonic levels predicted from Eq. A-21 are in excellent agreement with 
the data in Figs. D-l to D-5. The reason for the departure of the pre- 
diction from the measured data in Fig. D-6 may be. that the signal is too 
intense to be handled with the second order perturbation analysis. Some 
additional discussion of this point is given below.. 

Typical time waveforms observed during the high level experiment at 
5.33 kHz are shown in Fig. D-7. The progressive distortion of the wave with 
range is apparent. Referring to Fig. D— 6, we see that the second harmonic 
has grown 13 dB relative to the fundamental over the propagation distance from 
0.97 to 24.3 m. The asymmetry of the last three waveforms (sharp peaks, 
rounded troughs) is common and is caused by diffraction effects (see Ref. 1 
or Chapter 2). 
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FIGURE D-7 

MEASURED TIME WAVEFORMS FOR f = 5.33 kHz, 
EXPERIMENT AEM2. G STANDS FOR OSCILLOSCOPE GAIN. 
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An a priori estimate of the importance of nonlinear effects in the 
tone experiments may be found from the SFL chart shown in Fig. D-8 (see 
Chapter 2). The points shown were computed for the data given in Figs. 

D-l to D-6. The position of the points below the lower curve is consistent 
with our finding that nonlinear effects did not cause much extra attenua- 
tion of the findamental component. Notice, however, that the point for 
the high level experiment at 5.33 kHz falls on the lower curve, which 
marks the threshold of conditions for moderate nonlinear effects. It is 
therefore not surprising that the second order perturbation solution failed 
to provide a good fit for the second harmonic data in Fig. D-6. 

Beam patterns for the fundamental component were measured at a range 
of 18.2 m for both low and high level electrical input (Fig. D-9) . In each 
case the measured levels were normalized with respect to levels on the array 
axis. The solid curve is the pattern (computer by using linear theory) for 
a rectangular piston of dimensions equal to those of the array. Three 
comments may be made about the beam pattern measurements. First, 
the measurements were most difficult near the pattern nulls, where slight 
deviations in the beam caused by atmospheric inhomogeneities led to large 
changes in the received level. Second, except for a slight asymmetry, 
both the low and high level data agree well with the computed patterns 
for the major lobe. Third, the beam patterns measured at low and high 
levels differ at most by only 1 dB. Because the maximum difference in the 
low and high level data would be expected to be less than the extra attenuation 
suffered on axis for the high level signal, which has already been found 
to be miniscule, we conclude that nonlinearity had little effect on the beam 
pattern of the fundamental component. 
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0 * 



* LOW LEVEL - SPL = 126.7 dB 
O HIGH LEVEL - SPL = 144.7 dB 

COMPUTED FOR A RECTANGULAR 

PISTON (LINEAR THEORY) 


FIGURE 0-9 

MEASURED HIGH AND LOW LEVEL BEAM PATTERNS FOR 
AEM ARRAY AT 5.33 kHr. SOLID CURVE IS COMPUTED 
FOR A RECTANGULAR PISTON (LINEAR THEORY). 
MEASUREMENTS TAKEN AT A RANGE OF 18.2 m. 
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In summary, high and low level tests with tones were carried out 
with the AEM array. Propagation, waveform, and beam pattern measurements 
were made. Although considerable harmonic distortion developed, little 
extra attenuation of the fundamental component occurred in the farfield. 

A theoretical prediction for the second harmonic, where applicable, 
was found to be in good agreement with measured data. Beam patterns 
(at the fundamental frequency) at low and high levels were found to be 
practically the same. Both of these experimental finds are consistent with 
a priori estimates of the importance of nonlinearity based on use of the 
SFL chart. 
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